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Abstract 

We consider systems of interacting diffusions with local population regulation. Our 
main result shows that the total mass process of such a system is bounded above by the 
total mass process of a tree of excursions with appropriate drift and diffusion coefficients. 
As a corollary, this entails a sufficient, explicit condition for extinction of the total mass 
as time tends to infinity. On the way to our comparison result, we establish that systems 
of interacting diffusions with uniform migration between finitely many islands converge 
to a tree of excursions as the number of islands tends to infinity. In the special case of 
logistic branching, this leads to a duality between the tree of excursions and the solution 
of a McKean-Vlasov equation. 



1 Introduction 



In population dynamics and population genetics a prominent role is played by diffusion 
processes on [0, oo)'' or [0, 1]'' driven by stochastic differential equations (SDEs) of the form 



dXt{i) 



dt 



+ ^a^{Xt{i))dBt{i), ieG, (1) 



for t > where G is a finite or countable set and {Bt{i))t>Oi i & G, are independent standard 
Brownian motions and where {m{j, i))j^i(zQ is a stochastic matrix. We will refer to the 
solution {Xt)t>o of ([T]) as (G, m, /x, cr^)-process. In appropriate timescales and for suitable 
choices of /it and cr^, the component Xt(i) describes the (rescaled) population size on island 
j g G at time t > 0, or the relative frequency of a genetic type that is present on island 
« £ G at time t > Q. The linear interaction term on the right-hand side of ([T]) models a 
mass flow between the islands, which might be caused by migration of individuals or a flow 
of genes. Here we will use a picture from population dynamics. The coefficient <J^{x) then is 
the infinitesimal variance of the local population size given its current value x G [0,cx)). A 
classical case is Feller's branching diffusions where fT^(x) — const • x for x S [0, oo). Moreover 
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the drift term fj, describes the growth rate of a local population apart from immigration 
from other islands and emigration. A prototypical example is fJ-{x) = x{K — x), i.e. logistic 
growth. The migration matrix m(i, j)i^j^o is assumed to be a stochastic matrix. 

Here we address the question of the maximal effect of the migration matrix m{i, j)ij^G for 
fixed ^ and a^. We will show in Theorem [2] that the total mass process of (^4)00 is dominted 
by the virgin island model, which has been constructed in [19j . This dominating process does 
not depend on G or on the migration matrix. The intuition which leads to a comparison result 
is as follows. Consider a model with supercritical population-size independent branching and 
additional deaths due to competition within each island of G, i.e., /it is a concave function 
with 11(0) = and is a linear function. Now compare different distributions of individuals 
over space. If there is at most one individual on each island, then there are no deaths due 
to competition until the next birth or migration event. If, however, all individuals are on 
the same island, then there are deaths due to competition. The effect of the population-size 
independent branching is the same in both situations. We infer that more individuals survive 
if the distribution of individuals is more uniform over space. Now the migration dynamics 
which distributes mass uniformly over space would be uniform migration on G. As there is 
no uniform migration on an infinite set G, we approximate G with larger and larger finite 
subsets and consider uniform migration on the finite subsets. This intuition leads to the 
assertion that the total mass of the (G, m, cr^)-process is dominated by the total mass of 
an N -island model {X^)t>Q if we let — > 00 

N 

^A,(z)< hm ^Af(z), t>0. (2) 

zGG i=l 

For A e N := {1,2,...} the iV-island model {X^)t>Q is the solution of ^ with G := 
{1,...,A^} and m{i,j) :— i,j G G. We will refer to this A^-island model as {N,^,,a'^)- 
process. The initial configurations are assumed to converge suitably, as A^ 00, to a 
sequence {x{i))iiz^ whose sum is finite. 

Our first main result (Theorem [l]) establishes convergence of {XP)t>Q as A — > 00. The 
limiting process turns out to be a forest of (mass) excursions started in x{i), i g N. This 
forest of excursions has been constructed and analysed in [TH] and has been denoted as virgin 
island model. In this model, every migrant populates a new island, similar to the infinite- 
allele model in which every mutation produces a new allele. Jean Bertoin turned this idea 
into a tree of alleles and a partition into colonies [6] (see also Section 7 of Pardoux and 
Wakolbingcr (2010) for a connection with the virgin island model). In particular, Bertoin 
[7j shows in case of independent branching that the virgin island model is the diffusion 
approximation of a discrete mass branching process. In Section 2 we will review and discuss 
the virgin island model. 

Here is a brief heuristics how a forest of (mass) excursions emerges from the (A^, /i, cr^)- 
processes as A^ — > 00. For finite A^, every island continuously sends migrants to every 
other island. Every once a while (in fact, densely in time), a migrant is successful and 
found offspring which immediately find themselves interacting with the resident population 
on that island. For a bounded initial mass, the probability that two "large" populations 
that descended from two different founders coexist at the same time in one and the same 
island turns out to be negligible as A^ — 00. Thus, in the limit A^ — )■ 00, every emigrant 
populates a previously unpopulated island. The evolution of the population size on every 
freshly populated island is described by a random (mass) excursion. These mass excursions 
are born (densely in time) in a Poissonian manner on ever new islands with an intensity 
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proportional to the currently extant mass, and, once born, follow the SDE 



dYt - -Ytdt + fi{Yt)dt + ^aHYt)dBt . (3) 

Formally, this is described by means of the excursion measure Q associated with ([3| in the 
sense of Pitman and Yor (1982) (see also [19]). The intensity measure with which a path 
{zt)t>o spawns a "daughter" excursion born at time i > is Zt dtQ. The roots of the forest 
are random paths i g N, which are independent solutions of ([3| with yj'^*-*' = x{i), 

i g N. The virgin island model is then a countable familie V :— {(s,x)} of islands where 
island (s, x) is populated at time s > and carries mass Xt~s at time t > 0. 

One tree of excursions in the forest of excursions is illustrated in Figure [T] Note that Fig- 




Figure 1: Subtree of the Virgin Island Model. Only offspring islands with a certain excursion height are 
drawn. Note that infinitely many islands are colonized e.g. between times si and S2. 

ure[T]does not contain the whole tree. In fact islands are populated by emigrants densely in 
time but only finitely many excursions started by these emigrants reach a given strictly pos- 
itive height. Now a noteworthy observation is that the tree structure provides us with inde- 
pendence of disjoint subtrees. Putting this differently, the virgin island model is a branching 
process in discrete time in the sense of Jifina (1958) except that there are now infinitely many 
types, one for each excursion path. Due to this branching structure, the virgin island model 
is easier to study than the iV-island process. Several authors considered analogous of the vir- 
gin island model in the case of state-independent branching; see [Tl [51 1^ 171 [T^ IT51 IMl BT] 
for a selection of articles. 

To state our convergence result of Theorem [l] more formally, denote the population size 



3 



spectrum for the (TV, /i, cr^)-process as 



N 



2 = 1 

for every N Furthermore define the population size spectrum of the virgin island model 

as 

C*:- E ^x.-.^ i>0. (5) 
(s,x)ev 

Our convergence result assert that the population size spectra converge in distribution, i.e., 

^ C (6) 

in distribution as TV — >■ oo; see Theorem fTlfor the precise statement. 

Now we state the comparison result B more precisly. Recall that fi is the infinitesimal 
mean in a non-spatial situation. We assume /i to be subadditive. Then a population of 
size X that is separated into two islands experiences (in sum) a larger growth rate than 
a population of the same size that is concentrated on one island. Thus the virgin island 
model should offer in expectation a more prolific evolution of the total mass than a model 
([T]) with the same coefficients /i and tr^. The infinitesimal variance has an impact on a 
comparison in distribution. More precisely the stochastic order in ([2]) depends on whether 

is superadditive, additive or subadditive. In our prototype example of additive cr^, the 
stochastic order is the usual increasing order. Then we have that 



E 



[/(5]X,(z))] <e[/( ^ Xt-s)], t>0, (7) 



for every non-decreasing function /: [0,oo) -> [0, oo). If cr^ is superadditve or subadditive, 
then we will use a concave increasing order and a convex increasing order, respectively. In 
fact Theorem [2] is a comparison result not only on the one-dimensional, but on the finite 
dimensional distributions. Thereby results on the distribution of the total mass of the virgin 
island model have an immediate impact on the distribution of interacting diffusions with 
local population regulation. 

As a very special application of our comparison result, we obtain a sufficient condition 
for extinction for interacting diffusions with local population regulation. Here we speak of 
global extinction if the total mass '^i^Q Xt{i) converges to zero in distribution as t — > oo 
whenever J^iec -^oi'i) < c»- Theorem 2 of [19 shows that global extinction of the virgin 
island model is equivalent to 

;exp( / dx) dy<l. (8) 



As a consequence of Q, if condition ([s]) is satisfied, then the (G, m, /i, cr^)-process dies 
out globally no matter what the migration matrix m(-,-) is. Here are further implica- 
tions of our comparison result. Theorem 2 of [19] together with ([t]) implies an upper 
bound for the survival probability of {Xt)t>Q. Theorem 3 of fT9^ yields an upper bound 
for E[Jq°° SigG Xs{i)ds\ if the left-hand side of (|8]) is strictly smaller than one and an upper 
bound for the growth rate of '^^[2Zi<^G-^s(S)'is\ as t — >■ cx) if ([8| fails to hold. Moreover 
Theorem 4 of [TS] implies an upper bound for the growth rate of X^iec "'^t(*) as f — !• oo if ([s]) 
fails to hold. 
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The limit of the A^- island model is well-known if the random variables {X^{i): i < 
N} are exchangeable with common law z/, say. Here the situation is different due to the 
unbounded total mass. Emigrants still choose to migrate to island 1, say, with rate 
which is small in the limit iV — > oo. However, the total mass is of order N in the case of a 
nontrivial exchangeable initial configuration. Consequently there is a continuous immigration 
of migrants onto island 1 even in the limit N ^ oo. Moreover the dynamics of the A^-island 
model preserves exchangeability. Therefore the rate of immigration converges to the expected 
average over all islands as — )■ oo due to de Finetti's theorem. This leads to the following 
limiting process. The iV-island process {X^)t>o converges in distribution as — > cx) to 
(Mt(i): t > 0,1 £ N) where (Mt(i))f>o, i € N, are independent processes and solve the 
McKean-Vlasov equation 



see e.g. Theorem 1.4 in [40] for the case = 1 and /i bounded and globally Lipschitz 
continuous. Thus the components of the TV-island process become independent in the limit 
N oo. Kac (1957) described the first known instance of this phenomenon and called it 
propagation of chaos. After McKean (1967), propagation of chaos has then been intensively 
studied in the mathematical literature, see Sznitman's St. Flour lecture notes [ID]. Our 
convergence ^ is another instance of propagation of chaos. Here independence in the limiting 
object is independence of disjoint subtrees. 

In the prototype example of logistic branching, the McKean-Vlasov equation ^ and 
the virgin island model can be interpreted as forward and backward process. Corollary [2] 
shows for that case that the solution of the McKean-Vlasov equation and the total mass 
process are dual to each other. This duality derives from a self-duality of the A^-island 
model. This self-duality in turn derives from pathwise forward and backward processes in a 
graphical representation, see [3]. So in our prototype example, the virgin island model may 
be interpreted as backward process of the solution of the McKean-Vlasov equation. 

Here is a selection of models for which one could think of a similar comparison result. 
Mueller and Tribe (1994) investigate a one-dimensional SPDE analog of interacting Feller 
branching diffusions with logistic growth. Bolker and Pacala (1997) propose a branching 
random walk in which the individual mortality rate is increased by a weighted sum of the 
entire population. Etheridge (2004) studies two diffusion limits hereof. The "stepping stone 
version of the Bolker-Pacala model" is a system of interacting Feller branching diffusions 
with non-local logistic growth. The "superprocess version of the Bolker-Pacala model" is an 
analog of this in continuous space. Blath, Etheridge and Meredith (2007) study a two-type 
version hereof, which is a spatial extension of the classical Lotka-Volterra model. Fournier 
and Meleard (2004) generalize the model of Bolker and Pacala (1997) by allowing spatial de- 
pendence of all rates. A model in discrete time and discrete space is constructed in Birkner 
and Depperschmidt (2007). In that paper an individual has a Poisson number of offspring 
with mean depending on the current configuration and, once created, offspring take an inde- 
pendent random walk step from the location of their mother. 

2 The virgin island model 

The virgin island model without immigration has been introduced in |19j . Here we slightly 
generalize this model by adding independent immigration of mass. 

The virgin island model is an analog of ([T]) in which every emigrant populates a new 
island. Islands with positive mass at time zero evolve as the one-dimensional diffusion (yt)t>o 




(9) 
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solving (|3|. The following assumption guarantees existence and uniqueness of a strong [0, oo)- 
valued solution of equation ([s]), see e.g. Theorem IV. 3.1 in [5T]. 

Assumption A2.1. The set I is an interval of length \I\ which is either of the form [0, |/|] 
if \I\ < oo or of the form [0,oo) if \I\ = oo. The infinitesimal mean /i: / — >■ IR and the 
infinitesimal variance : I [0, oo) are locally Lipschitz continuous in I and satisfy ^{0) = 
= cr^(O) and if\I\ < oo, then ^(|/|) < = (T^(|/|). The function cr^(-) is strictly positive on 
(0, |/|) and the function /i(-) is globally upward Lipschitz continuous, that is, ^{x) — fJ-{y) < 
L^\x — y\ whenever x > y £ I where is a finite constant. Furthermore satisfies the 
growth condition cr^{y) < io-(y + y^) for all y £ I where L„ is a finite constant. 

Note that zero is a trap for (Yt)t>o, that is, Yt = implies Yt+s = for all s > 0. If we need 
the solution {Xt)t>o of ([T]) to be well-defined, then we additionally assume the migration 
matrix to be substochastic. 

Assumption A2.2. The set G is (at most) countable and the matrix {m{j,i))j^i^Q is non- 
negative and substochastic, i.e., m{j,k) > and X^ieG "^(-^i *) — ^ f'^^ all j,k £ G. 



Note that Assumptio nJA2. l| together with Assumption | A2 . 2] guarantees existence and unique- 
ness of a solution of ^ with values in {x £ l'^ : \x\ < oo}. This follows from Proposition 
2.1 and inequality (48) of 20J by letting (7^ 1 for i £ G and using monotone convergence. 

Mass emigrates from each island at rate one and colonizes new islands. A new population 
should evolve as the process (Yt)t>o and should start from a single individual which has mass 
zero due to the diffusion approximation. Thus we need the law of excursions of (Yt)t>o from 
the trap zero. For this, define the set of excursions from zero by 



C/:={xeC((-^,^),[0,oo)):Toe (0,cx)],xt=0 Vte (-oo, 0] U [Tq, cx))} 



(10) 



where Ty = Ty{x) ■— mi{t > 0: Xt = y} is the first hitting time of y G [0,oo). The set 
U is furnished with locally uniform convergence. For existence of the excursion measure Q 
and in order to apply the results of |19j, we need to assume additional properties of /x(-) 
and of o'^(-). For the motivation of these assumptions, we refer the reader to [in]. Assume 



/o 



dy < oo for some < e < |/|. Then the scale function S: I [0, oo) defined through 



"'^^^^-'dx), S{y) := r s{z)dz, z,y £ I. (11) 



is well-defined. 

Assumption A2.3. The functions /u(-) and cr^(-) satisfy 



dy < oo and 



\i\ 



''iy)siy) 



dy < oo 



(12) 



for some < e < |/|. 



Under Assumption A2.3 the process {Yt)t>o hits zero in finite time almost surely and the 



expected total emigration intensity of the virgin island model is finite, see Lemma 9.5 and 
Lemma 9.6 in 19j. Moreover the scale function S{-) is well-defined and satisfies S (0) — 1. 



A generic example which satisfies Assumption A2.3 is /i(y) — ciy^^ — C2?/"^, cr^(y) — c^y 



where ci,C2,C3 > 0, K2 > ki > 1 and K3 £ [1, 2). Assumption A2.3 is not met by cr^(?/) = y^. 
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Assuming A2.1 and A2.3 we now define the excursion measure Q as the law of {Yt)t>Q 
started in e > and rescaled with 5'(e) as e — J' 0. More formally, it is stated in Pitman and 
Yor (1982) (see [19] for a proof) that there exists a unique tr-finite measure Q onU such that 



lim^E- [i^(y)] = j F{x)Q{dx) (13) 

for every bounded continuous F: C([0, oo), [0,oo)) — > R for which there exists a. 5 > such 
that F{x) = whenever supj>o Xt < S. The reader might want to think of Q as describing 
the evolution of a population founded by a single individual. In the special case (j'^{y) = 2j3y, 
li{y) — cy with /3 > and c g R, the process (yt)t>o is Feller's branching diffusion whose 
law is infinitely divisible. Then the excursion measure coincides with the canonical measure. 

Next we introduce the state space of the virgin island model. The construction of the 
virgin island model is generation by generation. As 0-th generation, we denote the set of all 
islands which are non-empty at time zero or which are colonized by individuals immigrating 
into the system. All islands being colonized by emigrants from the 0-th generation are denoted 
as islands of the first generation. The second generation is colonized from the first generation 
and so on. For this generation-wise construction, we use a method to index islands which 
keeps track of which island has been colonized from which island. An island is identified with 
a triple which indicates its mother island, the time of its colonization and the population size 
on the island as a function of time. Let 

{(0,s,x) : s > 0,x e C([0,oo), [0,oo))} (14) 

be the set of all possible islands of the 0-th generation. For each n S N, define 

In := : hi-i e In-i, (s,x) e [0,oo) X C/} (15) 

which we will refer to as the set of all possible islands of the n-th generation. This notation 
should be read as follows. Island i„ = (i^n-iy SyX) ^ ha-s been colonized from island 
in-i G at time s and carries total mass Xt-s at time t > 0. Denote by I the union of 

all In over n S N>o. The virgin island model will have values in the set of subset of I. 

Having introduced the excursion measure, we now construct the virgin island model with 
constant immigration rate 6 £ [0,oo) and started in {xk)k£¥i C /. Let {{Yf'''^'')t>a: k S N} 
be independent solutions of ^ such that Y^'^'' = almost surely. Moreover let II" be a 
Poisson point process on [0, c») x U with intensity measure 

w.[n^{dt®dij)\^edt®Q{dii)). (16) 

The elements of the Poisson point process II" are the islands whose founders immigrated into 
the system. Next we construct all islands which are colonized from a given mother island. 
Let {LL : t e X} be a set of independent Poisson point processes on [0, oo) x U with intensity 
measure 

E [n(''^'>^) {dt (g) dij^)] = Xt-s dt (g) Q{d^p) t e I U {0}. (17) 

The elements of the Poisson point process n^^'**'''^) are the islands which descend from the 
mother island (t, s,x). All ingredients are assumed to be independent. Now the 0-th gener- 
ation is 

V(°) :={(0,O,r'='^'-): A:eN}u{(0,i,x): (i,x)en0}. (18) 

The (n -I- l)-st generation, n > 0, is the (random) set of all islands which have been colo- 
nized from islands of the n-th generation and the virgin island model V is the union of all 



7 



generations 

V("+i):-{(.„,t,x)eX„:.„eV("),(t,x)Gn'''} and V:= |J V^"). (19) 

nGN>o 

We call V the virgin island model with immigration rate and initial configuration (xfc)^^^. 

3 Main results 

We begin with convergence of the A^-island process. In this convergence, we allow the drift 
function /i^v and the diffusion function crjv to depend on N in order to include the case of 
weak immigration. For example, one could be interested in an TV-island model with logistic 
branching and weak immigration at rate ^ on each island. In that case, one would set 
^in{x) = ;^ + Jx{K — x) and (J%{x) := x ioi x £ I. The equation of the iV-island process 
now reads as 



1 ^ 



N ^ 



dt 



^al{X^{i))dB^{i) (20) 



where i — 1, . . . and where {Bt{i))t>o, i S N, are independent standard Brownian motions. 
The idea to include weak immigration into a convergence result is due to Dawson and Greven 
(2010) who independently obtain convergence of an iV-island model using different methods. 

Define jiN{x) :— ^,n{x) — /iiv(O) for x £ I. For the (A^, /i^r, cr^)-process to converge, we 
need assumptions on /xjv, cr^ and on the initial distribution. 

Assumption A3.1. Define I := [0,oo). The functions /iat,/^: / — R are locally Lipschitz 
continuous on I. The sequence {fJ-N)Ne'N converges pointwise to ^ as N ^ oo. In addition, 
N-hn{0) -> 6* e [0, cx)) as iV oo and < Nhn{0) < 26* for all TV € N. The diffusion 
functions af^ and are linear, that is, cr^(a;) = /S^x and o'^ix) — f3x for some constants 
PntP > 0. Furthermore /3jv converges to (3 as N oo. Assumptions \X2.1\ and \A2^ 
hold for fi. Moreover {iJ.N)Ne'N is uniformly upward Lipschitz continuous in zero, that is, 
fiN{x) — fJ-Niy) < L^\x — y\ for all x> y £ I , G N and some constant L^. 

Here is an example. If fiwix) = + Cix"^ — C2x'^^ and crfj{x) = x, then Assumption 
is satisfied if 1 < ki < K2 and C2 > 0. 

Assumption A3. 2. The random variables (Xo(i))igi\i and {X^ {i))i<:i\f are defined on the 
same probability space for each TV G N. There exists a random permutation 7rf , . . . , 7rj^ of 
{1, . . . , iV} for each TV G N such that 

N 

limEf^|Xo(z)-Xo^(7rf)|l =0. (21) 

i=l 

Furthermore the total mass of Xf){-) has finite expectation IE(^jgj^ Xo(i)) < 00. 

If (xi)igiv[ C / is a summable sequence, then Assumption A3.2| is satisfied for Xo(i) = Xi and 
aAT ■ 



A3.1 



X^{i) = Xi, i < TV e N. 

Next we introduce the topology for the weak convergence of the iV-island process. What 
will be relevant here is not any specific numbering of the islands but the statistics (or "spec- 
trum") of their population sizes, described by the sum of Dirac measures at each time point, 
that is, 



N 

^^X,»(^) ) (22) 
t<T 



i=l 
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where 5x is the Dirac measure on x. The state space of the measure- valued process (22) is 
the set M.f{I^ of finite measures on /. We equip the state space Mf{I) with the vague 
topology on / \ {0}. 

Now wc formulate the convergence of the (iV, /iat, CT^)-process defined in ( pO| ). 

Theorem 1. Suppose that {tJ-w) p^^-^ 



and (cr^ 



satisfy Assumption 



T e [0,cxd), we have that 

, N 



Hon 


A3.1 


ASA 


> Th 



and that the 



t<T 



(t,ji,77)ev 



as iV —> oo 



t<T 



(23) 



in distribution where V is the virgin island model with immigration rate 6 — limjv_ 
and initial configuration (Xo(i))igi\i- 

The proof is deferred to Section |4j 

Remark 1. For readability we rewrite the convergence in terms of test functions. The weak 
convergence ( 23 1 is equivalent to 

r / ^ 

^ (E/(^"«) 



hm E 



t<T 



= E 



F 



E 

(t,s,r()eV 



fivt- 



t<T 



(24) 



for every bounded continuous function F on C([0,r],R) and every continuous function 
/: / — ^ R with compact support in (0, |/|). Applications often require functions f with non- 
compact support. The following condition might be useful in that case. Let F be a continuous 
function on C([0,r],R) satisfying the Lipschitz condition 



\F{ri) 



<^fEI% 



V7y,7^€C([0,T],R) 



(25) 



for some < ii < • • • < i„ < T. In addition let f : / — > R be a continuous function satisfying 
|/(a:)| < Ljx for all x £ I. Then following the arguments in the proof of Lemma 22 below, 
one can show that (24) holds with F and f replaced by F and f, respectively. o 

The assumptions of Theorem[T]are satisfied for branching diffusions with local population 
regulation. A prominent example is the A^-island model with logistic drift fi{x) = jx{K — x) 
and with cr'^{x) = 2/3x for x E [0, oo) and some constants 7,X, /3 > 0. More generally, 
Theorem [l] can be applied if /i(x) = 7a; — c{x) and cr^(a;) — 2/3x for x G [0,oo) where c(-) 
is a concave function with c (0) E R. We believe that Theorem [l] also holds for non-linear 
infinitesimal variances such as cr^(a;) = a::(l — a;) in case of the Wright-Fisher diffusion. Our 
proof requires linearity only for one argument which is the step from equation (144) to 
equation (145). 

In case of logistic branching, we obtain a noteworthy duality of the total mass process 

Vf.^ J2 Xt-s, t>0. (26) 
('.,s,x)ev 

of the virgin island model with the mean field model {Mt)t>o defined in ([9]). By Theorem 3 
of [20], systems of interacting Feller branching diffusions with logistic drift satisfy a duality 
relation which for the {N, ^y{K — y), 2/32/)-process reads as 



'(i)y 



Vx,y,i>0 



(27) 
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where the notation Wj'"^^ refers to the initial configuration — {y,0, . . . ,0) and JS^- refers 
to = {x,. . . ,x). This duality is established in [2(7 via a generator calculation, in Swart 
(2006) via dualities between Lloyd-Sudbury particle models and in j3] by following ancestral 
lineages of forward and backward processes in a graphical representation. Now let — >■ cxd 
in (27). Then the left-hand side converges to the Laplace transform of the total mass process 
of the virgin island model (without immigration) and the right-hand side converges to the 
Laplace transform of the mean field model This proves the following corollary. 

Corollary 2. Let {Vt)t>o be the total mass process of the virgin island model without im- 
migration starting on only one island. Furthermore let {Mt)t>o be the solution of both 
with coefficients /Lt(y) = jy{K — y) and cr^(y) = 2(3y for y £ [0, oo) where j,K, /S > 0. Then 

ES/e"?^^' =E^e-p*^'^ Vx,y,i>0 (28) 
where and refer to Vo = y and Mq = x, respectively. 

Together with known results on the mean field model, this corollary leads to a computable 
expression for the extinction probability of the virgin island model. 

Corollary 3. Let {Vt)t>a be as in Corollary^ Then Vt converges to a random variable Vaa 
in distribution as t ^ oo. Lf 

exp f K^x —x^ j • exp(— .x) dx < 1, (29) 



then P(Voo = 0) = 1. // condition (29 I fails to hold, then 



oo 

" yx - 



]py{V^ = 0) = 1 - pyiV^ = oo) = / e"?^" Tpidx) e (0, 1) (30) 

Jo 

where refers to Vq ~ y G {0, oo). The parameter p € (0, cxi) is the unique solution of 

y^K - y) e^p {^^^^y ~ ^y^^y ^ (31) 
and the probability distribution Fp is defined by 



Tpidx) = exp / '-^ '—^ '-dz dx (32) 

px \Jk pz J 

on (0,cx)) where Cp is a normalizing constant. 

Proof. Theorem 2 of JH] shows convergence in distribution of Vt to Voo as i — >■ oo and 
V{Voo = 0) = 1 if ([29]) holds. If ([29| fails to hold, then Corollary[2]together with convergence 
of Vt implies convergence in distribution of Mt to a variable Moo as i — >■ cx). The distribution 
of Moo is necessarily an invariant distribution of the mean field model (|9| and is nontrivial. 
Lemma 5.1 of shows that there is exactly one nontrivial invariant distribution for ([£]) 
and this distribution is given by ( 32 ) . □ 

The second main result is a comparison of systems of locally regulated diffusions with 
the virgin island model. For its formulation, we introduce three stochastic orders which are 
inspired by Cox et al (1996). Let Z = {Zt)t>o and Z = {Zt)t>o be two stochastic processes 
with state space /. We say that Z is dominated by Z with respect to a set F of test functions 
on path space if 

Z <F E[i^(Z)] < E[F(Z)] Vi^eF. (33) 
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The first order is 'tlie usual stocliastic order' <st in which Z is dominated by Z if there is a 
coupling of Z and Z in which Zt is dominated by Zt for all t > almost surely. Assuming 
path continuity, an equivalent condition is as follows. Denote the set of non-decreasing test 
functions of rt e N>o arguments by 

J"|'2 := -7^|±(/) := |/: /" ^ ]R| / is non-decreasing, / is bounded or / > o| (34) 

for a set S C [0,oo). Furthermore let F+j- be the set of non-decreasing functions which 
depend on finitely many time-space points 

F+± :=F+±(G,5) 

:= {f: C([0,c») x G, 5) ^ R|37i G Nq 3 (ti, zi), . . . , (t„, i„) £ [0,oo) x G (35) 
3/e J-|"2(^) such that F{r])^ f{r,t,{ii),...,vtjin))}. 

If there is no space component, then we simply write F_|_-|-(S'). In this notation, Z <st Z is 
equivalent to Z <f+± Z, see Subsection 4.B.1 in Shaked and Shanthikumar (1994). 

We will use two more stochastic orders. In the literature, the set of non-decreasing, 
convex functions is often used. Here an adequate set is the collection of non-decreasing 
functions whose second order partial derivatives are nonnegative. As we do not want to 
assume smoothness, we slightly weaken the latter assumption. We say for 1 < i,j < n that 
a function /: R" — ^ IR is (i, j)-convex if 

/(z + h^e, + h^ej) - f{z + /lie,) - f{z + h2ej} + f{z) > V z G R", hi,h2> 0. (36) 

Note that if / is smooth, then this is equivalent to g^. / > 0. In addition note that / 
is (i,i)-convcx if and only if / is convex in the i-th component. Moreover we say that / is 
(i, j)-concave if — / is (z, j)-convex. A function is called directionally convex (e.g. Shaked 
and Shanthikumar 1990) if it is (j,j)-convex for all 1 < i,j < n. Such functions are also 
referred to as L-superadditive functions (e.g. Rueschendorf 1983). Define the set of increasing, 
directionally convex functions as 

J"|"2 := { / e ^+1 ■ I is (i,i)-convex for aU 1 < j < n| (37) 

and similarly J--\ with '(i, j')-convex' replaced by '(z, j)-concave'. Furthermore define F_|_+ 

and F_| as in (351 with replaced by J-"^"^ and respectively. Now we have intro- 

duced three stochastic orders <f+±, ^f+_ and <f++- Note that contains all mixed 

monomials and that F^^l contains all functions 1 — exp (— X^iLi ^i^i) with Ai, . . . , A„ > 0. 



Theorem 2. Assume A2.1, A2.2 and A2.3. If ^ is concave and if cr^ is superadditive, then 
we have that 

(E^^«) <F.-ao,oo), (38) 

///Li is concave and is subadditive, then inequality (38) holds with F_| replaced by F_|__|_. 

If H is subadditive and is additive, then inequality (38) holds with F_| replaced by F-|_±. 

The proof is deferred to Section [5j Comparisons of diffusions at fixed time points are well- 
known. Theorem [2] provides an inequality for the whole process. The techniques we develope 
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for this in Subsection |5.1| might allow to generalize the comparison results of Bergenthum 
and Rueschendorf (2007) on semimartingales. 

The assumption of /x being subadditive is natural in the following sense. Let us assume 
that letting two 1-island processes with initial masses x and y, respectively, evolve indepen- 
dently is better in expectation for the total mass than letting one 1-island process with initial 
mass X -\- y evolve. This assumption implies that 



u[x + y) = hm 

t-i-O 



< lim 

t t^o t 



x-y 



ii{x) + M(y) (39) 



for all X, y, a; + 2/ G / and thus subadditivity of the infinitesimal mean /i. If is not additive, 



then we need the stronger assumption of /i being concave for Lemma 30 

From Theorem [2] and a global extinction result for the virgin island model, we obtain 
a condition for global extinction of systems of locally regulated diffusions. According to 
Theorem 2 of , the total mass of the virgin island model converges in distribution to zero 
as t — )■ oo if and only if condition (40 1 below is satisfied. Together with Theorem 2, this 
proves the following corollary. 



A2.1 


A2.2 


and 


A2.3 



Corollary 4. Ass'^ 
or that /i is concave and is either superadditive or subadditive. Then 



o\y)l2 



exp 



a^{x)/2 



dx] dy < 1, 



implies global extinction of the solution {Xt)t>o of ([T]), that is, X]ieG"^*(^) 
whenever X^ieG '''^o(i) < almost surely. 



additive 



(40) 



as t ^ oo 



In case of logistic branching {iJ.{y) = jy{K — y), a'^{y) ~ "^fSy), condition (40) simplifies to 
condition (29 1. For a second example, let {Xt)t>o be a stepping stone model with selection 
and mutation, e.g. / = [0, 1], /i(y) = sy{l — y) — uy and <j'^{y) — 2y(l — y) where s, u > 0. It 
is clear that (Xt)t>o dies out globally if s < u. Corollary |4] now improves this so that 



/ (l-y)"exp(sy)dy< 1 

JO 



(41) 



implies global extinction. 



4 Convergence to the virgin island model 
4.1 Outline 

First we outline the intuition behind the proof. The virgin island process is a tree of excur- 
sions whereas the iV-island process has no tree structure. It happens in the latter process 
that different emigrants colonize the same island. In addition, the A^-island process is not 
loop-free. An individual could migrate from island 1 to island 3 and then back to island 1. 
That these two effects vanish in the limit as the number of islands tends to infinity will be 
established in two separate steps. 

The first step ensures that the limit of the iV-island process as iV — ?> oo is loop-free. 
For this purpose, we decompose the A^-island process according to the number of migration 
steps. Throughout the paper, we say that an individual has migration level k € Ng at 
time t e [0, oo) if its ancestral lineage contains exactly k migration steps. For example. 
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an individual starting on island 1 at time 0, moving to island 3 and then back to island 
1 has migration level 2. Let iV G N. We define a system {{X^'''{i))t>o: k G No,i < N} 
of diffusions such that ' (i) consists of the mass on island i at time t with migration 



level k. Recall fj,N{x) — I^n{x) — MAr(O) for all x £ I. Define Xq ' (i) — lk=ciXQ (i) for all 
i e {1,. . . ,iV} and fc e Nq. Let {{X^'''{i), B'^{i))^^^: i < N,k e Nq} be a solution of 



dxn^) = 



f.^(j2x^-"^{{))dt 



m>0 



(42) 



Xt W 



m>0 Xt 



N.i 



'N 



m>Q 



for all t G [0, oo), i — 1, . . . , N and k G Nq where X^'^^ := for i > and i < N and where 
{{B^{i))t>o: i,k G Nq} is a system of independent standard Brownian motions. Here we 
implicitly used the continuous extension of -^^rpyfiN{x + y) and of -^^r^crl^ix + y) as functions 
of {x,y) e [0, oo)2 \ {(0,0)} into the point (0,0), where iV e N. Any weak solution of (|42]) 



will be denoted as a (iV, /xat, cT^)-process with migration levels. See Lemma [H] for existence 
of a weak solution of (42). 



Lemma 24 below indicates that the individuals with migration level fc at a fixed time are 
concentrated on essentially finitely many islands in the limit N oo. A later migration 
event will not hit these essentially finitely many islands because hitting a fixed island has 
probability . Therefore we expect that all individuals on an island have the same migration 
level. Inserting this into (42) suggests to consider the solution {{Z^''^ {i))t>o : i < N,k G Nq} 
of 



dzn^) 



1 ^ 



t u) 



(i) + lfc=oMjv(0) + Ajv(^t ' (0) dt 



(43) 



r%{Z^''i^))dBn^), 



ryN,k^-\ vN.k,-\ 



where Z^' ^ := for all t > and i < N. We will refer to this solution as the loop- 
free (A^, /iAT, cr^)-process or as loop-free A^-island model. Lemma 26 below establishes the 



assertion that the distance between the (A", /zat, (T^)-process with migration levels and the 
loop- free (A^, /iat, cr^)-process converges to zero in a suitable sense as A — > cx). It turns 
out that some difhculties arise from the different forms of the diffusion coefficients in the 
(A', /ijv, (T^)-process with migration levels and in the loop-free (A^, ^jy, o-^)-process. As we 



could not resolve these difficulties, we additionally assume for Lemma 26 that aj^ is linear. 



Then we have that xcj'j^{y)/y — cr^(x) for x — X^'^{i) and y — Yl°^=oX-. 



diffusion coefficients in (^42h and in (43) are similar. Our proof of Lemma 26 



rN,m 
'■t 



(i) and the 
is a moment 

estimate in the spirit of Yamada and Watanabe (1971). 

In Subsection |4.3| we show that two emigrants colonize different islands in the limit A^ — !• 
oo. Let us rephrase this more formally. Recall that {{B^{i))t>Q: i < N} is independent of 
{{Bi{i))t>o: I < k,i < N}. Thus, conditioned on {E^i '''"^(j) : t > 0}, the loop-free 
(A^, /iAT, (T^)-process with migration level fc > 1 evolves as the solution of 



dYZH^) 



A^ 



dt-Y, 



!;';H^)dt + Iin{y/'-H^)) dt + ^fjAY^^dBU^), (44) 
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where Y^^/{i) = 0, s = 0, and CtvW := T,f=i Z^'''~\j)- Note that {Yt^;Hi))t>s, i < N, are 
independent and identically distributed. Now we are interested in the total mass li^^'*"] as 
iV ^ oo. As Yf^'^{i) = and as the immigration rate on island i tends to zero, the process 

{Yf^g'^ {i))t>s converges to zero as iV — )• oo for every i g N. However, as mass of order O(^) 
immigrates on a fixed island, the probability that the excursion started by these immigrants 



reaches a certain level 5 > 0, say, is of order O(^) as the convergence in (131 indicates. Now 



as there are N independent trials, the Poisson limit theorem should imply that 



N 



as N oo 



(45) 



/ s<t<T \J / s<t<T 

where 11 is a Poisson point process with intensity measure 



N 



lim (N{u)du Q{dri). 



(46) 



We will prove (451 in Lemma 22 by reversing time. 

For convergence of the loop-free (N, tT^)-process, we do not need to assume linearity 
of the diffusion function. Here we may replace Assumption |A3.1| with the following weaker 
assumption. 



Assumption A4.1. The functions /iAr,/i: / 



7|^,cr^; / — > [0,oo) are locally Lips- 



chitz continuous on I . The sequence (i^n, <^m)n<£'N converges pointwise to {fi, a ) as iV — > oo. 
In addition, N-fiN{0) ^ 9 € [0, oo) as N ^ oo and 7V^Ar(0) < 20 for all N elN. The func- 
tions fiN anda% s atisfy /xat (0) > = cr^(O) and if \I\ < oo, then /i^vd^l) < = (t^(|/|). As- 
l A2.3 hold for ^ and . Moreover {^N)NeK is uniformly upward Lipschitz 



A2.1 



and 



sumption 

continuous in zero, that is fj,f^{x)—iJ,N{y) < Lfj_\x — y\ for all x > y £ I , N € ¥! and some con- 
stant L^. The sequence ((T^)ArgN satisfies the uniform growth condition crff{y) < L„(y -\-y'^) 
for all y I, N £¥! where L„ is a finite constant. 



A4.1 



implies Assumption A3.1 



Note that if af^ is linear, then Assumption 

Some steps of our proof are based on second moment estimates and require the following 
assumption of uniformly finite second moments of the initial distribution. This assumption 
is then relaxed in further steps. 

Assumption A4.2. The initial distribution satisfies that 



N 

supEffV<(*) 



=1 



< oo. 



(47) 



4.2 Preliminaries 

In this subsection we establish preliminary results such as moment estimates and existence 
of the processes. The quick reader might want to skip this subsection. We begin with weak 
existence of the A^-island process with migration levels. 



Lemma 5. Assume A4-.1 The (N, ^^,a'j^)-process with migration levels exists in the weak 
sense, that is, equation (42 I has a weak solution for every iV € N. 



Proof. As the proof is fairly standard, we only give an outline. Approximate (42) with 



stochastic differential equations for which weak solutions exist. For example, approximate 
fj.N and a% with bounded continuous functions /i7v,i 



and crfj n, respectively. Consider the 
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solution {X^'^'^)t>Q of (42 1 with and cr^ replaced by ^N,n and cr^„, respectively, and 
which only depends on the migration levels k < n. Then this solution has a weak solution 
according to Theorem V.23.5 and Theorem V.20.1 of [34] as the coefficients are bounded 
and continuous and the stochastic differential equation is finite-dimensional. Show that the 
formal generator hereof converges to the formal generator associated with (42). In addition 
establish tightness of Then there exists a converging subsequence and its limit 

solves the martingale problem associated with (42), see Lemma 4.5.1 in [15]. From this 
solution of the martingale problem, construct a weak solution of (42) as in Theorem V.20.1 
of 134 . □ 



Next we prove that the TV-island model with migration levels is indeed a decomposition 

Fix N elN and let { (xf '''(i), (i))^^^ : kelNo,i<N} be a 



of the iV- island model pO). 
Lemma 6. Assume 



A4.I 



solution of (42). Then {{X^ {i))t>o: i < defined through 

t>0,i<N, 



fc>0 



(48) 



is the unique solution of the N -island model ( 20 ) corresponding to standard Brownian motions 
defined through 



dBtH) 



'(■0>o X! 

fc>0 



vN,k/ -N 



dB^^{i) + lt>r^dB]{i) i>0 (49) 



for every i < N where tn := inf{t > : Y^iLi I]m>o ^(^'"(O = 0}- 

Proof. The process {Bt{i))t>o is a continuous martingale with quadratic variation process 
[B{i), B{j)]t = 5ijt. Therefore Levy's characterization (e.g. Theorem IV. 33.1 in [34|) implies 
that ( 49 ) defines a standard Brownian motion. Moreover it follows from summing ( 42 ) over 
k g No that (X(^(i))^^p solves (20). Pathwise uniqueness of (20) has been established in 
Proposition 2.1 of [20]. ~ □ 



For the following lemmas, let { (X^ ' {i),B^{i))^^^ : fc € Nq, z < iV} be a solution of ^ 
and let {(^Z^''' (i)) ^^^i k e No,i < be the solution of (43). Define the stopping times 
£ [0,00) through 



N 



(50) 



i=l m>0 



for every KG [0,oo) and every TV S N. 



Lemma 7. Assume 
that 



lAjl] Then, for every T < 00, there exists a constant Ct < 00 such 



supE^ 

t<T 



N 00 



1=1 k=0 
N 



N 



(51) 



for every initial configuration Xq — x G / and every N € IN 
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Proof. Fix N e and T < oo. By Assumption A4.1 we have that fJ.N{x) < L^x + ^ for 
all X £ I. According to Lemma |6j {(X/^(z))t>o : * < defined through (48) is a solution 
of (20). Sum (20 1 over i < N, stop at time and take expectations to obtain that 



E 



2eds 



(52) 



for every t < T and K G [0, oo). Note that the right-hand side is finite. Now Gronwall's 
inequality implies that 



E 



N -, N 

4=1 J j=l 



(53) 



for alH < T and K E [0, oo). Letting K oo Fatou's lemma yields that 



supE 

t<T 



N 



'-i=l 



E (0 = sup E ^ hm inf (z) 



N 



1=1 

< supliminf E 

t<T K~¥oo 



N 



N 



(54) 



This proves inequality (51 1 with Ct ■= T ejq){LT). The inequality for the loop-free A^-island 
process follows similarly. □ 



Lemma 8. AssuTne \A3.^ and \A4.1\ Then we have that 



N 

E sup supE (^X^'\^)+^Z^'''{t 



i) < oo 



(55) 



for every T e [0, oo). 



Proof. We prove inequality ( 55 ) for the solution of ( 42 1 . The estimate for the solution of ( 43 ) 
is analogous. 

Recall flN{x) — /iAr(a;) — ^Af(O) for all x £ I. Apply Ito's formula to Y^^=i Xt^''^{i), take 
expectations, estimate //at (a;) < Lf^x for all x £ I and take suprema to obtain that 



N 



N 



sup supE^xf''=(i) < lk=o sup E^X^( 



WeN t<T 



i=l 



[ 2eik=o + sup supE V (x^^'"\i) + L^Xf ''^(i)) dt 

Jo N£¥! s<t ~! \ ) 



(56) 



for all T > and all fc € Np. Note that the right-hand side is finite due to Lemma [7] and 
Assumption |A3.2[ Summing over k < K £Y{ and applying Gronwall's inequality implies 
that 

K N M 

y sup supEV ''^(z) < ( sup ISyX^'ii) + 2eT) ■ e(i+^.)^ (57) 



for every if S N. Letting A' — > oo proves (55 1. 



□ 
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Lemma 9. Assume\JJ^ and\XJ^ Th 



en we nave 



that 



sup E 

JVSN 



N 



< oo 



(58) 



is finite for every T < oo. The analogous assertion holds for the loop-free {N, fi^^aj^)- 
process. 

Proof. Recall from (50) and fix iV e N for the moment. Acco rding to Lemma [6j 



20 



N.i 



is an iV-island model. Recall from Assumption 



A4.1 



that ijln{x) < L^x- 



=: flivix) for all x £ I. Thus Lemma 3.3 of [2^ imphes that the (A^, //at, cr^)-process is 
dominated by the (A^, fiN, cr^)-process (A/^)t>o. Using Ito's formula and <J%{x) < Lc{x+x'^) 
for all x € /, we get that 



E 



^ i=l 
t 



E 



< E 



2i,,E 



N 



1=1 



i=l 
2 



ds (59) 



1=1 

t r N 

(2L^ + L,)E 







+ (46* + L„)Ct 



for every t < T, every AT S N and some constant G [l,oo). We used Lemma [7| for the 
last inequality. Note that the right-hand side is finite. Applying Doob's submartingale 
inequality (e.g. Theorem n.70.2 in [33 ) to the submartingale X]iLi^t^(*)i using Fatou's 
lemma and applying Gronwall's inequality to (59), we conclude that 



E 



N 



t<T ^„ 

— 1=1 m>0 



AT 

<4supE[(5]Af(z) 

N 



< 



4 sup E [ f V lim inf Ai^ (i) 

— i—1 

N 

4supliminfE[fy Ai);~„(i; 

— i—1 

N 

+ L.)CtT + e(^Ao^( 



< 4 



i=l 



The right-hand side is bounded uniformly in A' € N due to Assumption |A4. 2} The proof in 
the case of the loop-free A^-island model is analogous. □ 



Recall from (50). Next we show that stopping at the time fj^ has no impact in the 
limit A' — !■ oo. 

Lemma 10. Assume 



A^.l and A4-.2 Then any solution of (42 1 satisfies that 
lim sup sup El sup^^ A(^''"(i)l^]v<rp j = 

^ ~ 1=1 m>0 ' 



(60) 



for every T € [0,oo). The analogous assertion holds for the loop-free {N, ^N,afj) -process. 
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Proof. Rewriting {f^ < T} = {supt<2- X^i^i Sm>o "'^t^''"(*) — ^^-^ assertion follows 
from the Markov inequality and from the second moment estimate of Lemma |9] □ 

Next we prove some preliminary results for the solution (yf^'^)f>o of ( |44[ ). 

;e A4-.1 Let (^n: [0,oo) — > N-I be locally square Lehesgue integrable 



Lemma 11. Ass'i 
function. Then we have that 



sup Y^j'^j 

s<t<T 



< Ci 



N 



dr 



(61) 



for allx€l,0<s<T, N(z¥! and some constant Ct < oo which does not depend on x, 
N or on C.pf. 

Proof. The proof is similar to the proof of Lemma [9j so we omit it. □ 

A4.1\ Let (i^ : [0, oo) — ?> N-I be locally square Lebesgue integrable 
i < N, be independent solutions of (44 1 for every 



Lemma 12. Assume _ 
function. Furthermore let (^^^'^(i)) 
iV e N. Then we have that 



sup 

,s<t<T 



N 



1=1 



< Ct 



C,N{r) + — dr 



(62) 



for all < s < T , N ¥! and some constant Ct < oo which does not depend on N or C,iq . 
Proof. The proof is similar to the proof of Lemma [9j so we omit it. □ 

Lemma 13. Assume 



A4.I 



and fixT ^ [0, 00). Let {Y^g^)t>s and (i^(^''')t>s be two solutions 
of (44) with respect to the same Brownian motion such that Y^^'^ = x and Yj^^''' = y. If 
Cn,Cn- [0,r] — N-I are square Lebesgue integrable, then for some constant Ct < 00 being 
independent of (^N,^N,x,y and N £ ¥! we have that 



E 



N..C 



Y 



N,C\ 



< Ct 



\CN{r) - CN{r)\dr+ \x - y\ 



(63) 



for all N e'N, for alio < s <t <T and all x,y e I. 



Proof. As in Theorem 1 of Yamada and Watanabe (1971) yt2j, an approximation of x \x\ 
with C^-functions results in 



dY, 



NX 



Y 



N.C\ 



sgn {Yt 



N..C 



Y 



Y 



Taking expectations, the upward Lipschitz continuity of /In implies that 



E 



NX 



Y 



NX\ 



<\x-y\ + j^ I \CN{r) - CN{r)\dr + L^ 



]&\Yl'/ -Yl'/ldr 



(64) 



(65) 



for all t>s. The right-hand side is finite due to Lemma 11 Therefore, Gronwall's inequality 
implies ^ with Ct = e^'^'^ . □ 
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Now we study the solution {Y/^g''^)t>s of (44) in which = c E N-I. Let a denote a 
fixed point in (0, |/|). Recall the scale function 5* from (11). Define the scale function Sn of 
iK^s'')t>s through 



sn{z) 



exp I 



c/N 



■ dx 



<^%{x)/'2 



dx), SNiy) 



for all y,z E I. The next two lemmas involve the speed measures 

2 . , . . 2 



m{dy) 



dy, mM{dy) 



<^%iy)sN{y) 



dy 



spf{z)dz (66) 



(67) 



as measures on / for every S N. Note that under Assumption A4.1 mAr(-) converges as 
measure on (0, |/|) vaguely to m(-). 



Lemma 14. Assume 
integrable, then 



A4.I 



and 



A2.3 If C,N,C,N'- [0,00) — > N-I are locally square Lebesgue 

mNidx)=0 (68) 



limlimsup / M'^fiY,^,''^) - ^'f{Y,^;^) 



for all s < t < 00 and all globally Lipschitz continuous functions f : [0, 00) — >■ R. 

Proof. Fix s <t and let C £ [0,oo) be the global Lipschitz constant of /. The first moment 
estimate of Lemma [13] provides us with the inequality 



lim sup 

N^oo Jo 



lE^f m^idx) 



(69) 



no po 

< lim sup / CxmM{dx)~ I Cxm{dx) 
TV— >oo Jo Jo 

for every S > 0. The last step follows from the vague convergence of mj\j to m together with 
dominated convergence and Assumption A2.3 The right-hand side of (69) converges to zero 
as (5 by the dominated convergence theorem together with Assumption |A2.3 □ 

and 



Lemma 15. Assume 
then 



A4.I 



A2.3. Let rriN be the speed measure of {Y/^g''^)t>s- lfc>0, 

(70) 



TV 



ruNmS)) 



1 



as iV —)■ 00 



for every S > 0. 

Proof. Fix S > and c > 0. Integration by parts yields that 



N 



mw((0,<5)) 



c/N 



<(2/)/2 



exp 



( - 



exp 



c/N 



^ ~x + p-wix) 
<(2;)/2 



dx] dy 



a "|/(2^)/2 



dx 



^ —X + jj,N{x) 

crI/(a;)/2 



dx 



1 s 



(71) 



^v + P-Niy) 



t^^(y)/2 



exp 



( - 



c/N 



a '^Ar(2^)/2 



dx 



^ —X + jj,N{x) 

^Nix)/2 



dx ) dy 
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for every TV €E N. As aff is Lipschitz continuous in [0, a], g-a ^^^2 '^^ ~ —00. Letting 
A'^ — >■ cx) in ( 71 ) and applying the dominated convergence theorem shows that 



exp 



exp 



, no 



— X + ^{x) 



dx I — 



dx 



^'(y)/2 



exp 



— 2/^^ /o ) d-y 



exp 



/o ^^(^)/2 
which is equal to one. 

We recall the following lemma from [TH], see Lemma 9.8 there. 





^ _^ s 

dx 



(72) 



□ 



Lemma 16. Assume A2.1 and A2.3 Let Q he the excursion measure defined through (13|. 
Then 



Xt dt Q{dx) 



\i\ 



y 



^'(y)/2 



exp 



^ ~x + iJ.{x) 
a^{x)/2 



dx ] dy < 00. 



(73) 



The last result of this subsection is a variation of the second moment estimate of LemmaO 
Define the stopping times G [0, 00) through 



N 

r#:=inf{t>0: ^Af(.)>if} 



(74) 



i=l 



for every KG [0,oo) and every iV S N. 
Lemma 17. Assume\IJ^\Ml\and\Il^ Th 

N 



sup 



en we have that 

2n 



t<T ^ 



< 00 



(75) 



is finite for every T < 00 and AT € N. 

Proof. Lemma 3.3 in [20] shows that, on the event {t^ > t}, (i) is stochastically bounded 



above by q' 



N,K+Niin(0) 



By Assumption 



A4.1 



Together with the second moment estimate ot 



we have that NfiN{0) < 20 for aU A^ e N. 
jcmma 



11 this implies that 



N 

1=1 



supYto 

t<T 



2n 



<CTJ2(M[X^i^)]+M[[x^{{)) 

i=l 

< Ct f sup EIAq*^! + sup E 



MSN 



K + 20 {K + 2ef 
N A2 

T{K + 2e)+T{K + 2ef 



(76) 



for every A' € N and some constant Ct < 00. The right-hand side is finite due to Assump- 
tions and [l42l □ 
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4.3 Poisson limit of independent diffusions with vanishing immigra- 
tion 

In this subsection, we prove (45) which is the central step in the proof of Theorem [I] Our 
proof is based on reversing time in the stationary process. For the time reversal, we consider 
the following stationary situation. Excursions from zero of the process {Yt)t>o start at 
times given by the points of an homogeneous Poisson point process on R with rate 1. This 
process of immigrating excursions is invariant for the dynamics of (Yt)t>o restricted to non- 
extinction. Now the time reversal of an excursion is again governed by the excursion measure, 



see Lemma 18 As a consequence reversing time in the process of immigrating excursions 
does not change the distribution. 

Let us retell the story more formally. Consider a Poisson point process 11 on (— oo, oo) x 
U with intensity measure ds Q. Then J2(s j))en^(vt-s)t>o process of immigrating 

excursions. Note that at a fixed time t, {jyt-s : (s, t]) € 11} is a Poisson point process on 
(0, cxd) with intensity measure 

/oo /"O 
QiVt-s e dy) ds = / Q(77_, e dy) ds = m{dy) (77) 
oo J —oo 



where m is the speed measure defined in (67). Here we used that r]t — for t e (— c)o,0]. 



The relation ( 77 ) between the speed measure and the excursion measure has been established 



in Lemma 9.8 of by exploiting a well-known explicit formula for /q f{Ys)ds. It is 
also well-known (e.g. equation (15.5.34) in [53]) that the speed measure m is an invariant 



measure for the sub-Markov semigroup 'E' f(Yt)lYt>o- This fact can also be seen from (77) 
by noting that Q(?7s G dy) is an entrance measure for this sub-Markov semigroup. Thus the 
process of immigrating excursions is indeed invariant for the dynamics of (Yt)t>o restricted 
to non-extinction. 

Now we show that reversing time in the process of immigrating excursions does not 
change the distribution of the process. The process {Yt)t>o restricted to non-extinction is 
time-reversible when started in the invariant measure m, that is, 

J M-F{iYt),^^)lY^yom{dx) = J M-F{iYT-t)t<T)'^Y^>om{dx) (78) 

for every T € [0, oo) and every non- negative measurable function on C([0,r]). If the speed 



measure m is replaced by the left-hand side of ( 77 ) , then ( 78 ) can be extended to allow for 
extinction. To show this we first formulate the Markov property of the excursion measure. 
Definition (13) of Q as rescaled law of {Yt)t>o together with the Markov property of (Yf)t>o 
implies that 

Q<t<T 

(79) 

= J F{{vt)o<,<T)^'^F{{Yt),^,)Q{dv) 

for all F,F: C([0,oo), [0,oo)) satisfying F{0) = = #(0) and every T e [0,oo). 
Lemma 18. Assume [A2l\ and [AE3{ Then 



P{irit-s)tei,)dsQidv)^ F{ivT-t-s)teTddsQidii) (80) 

) J J —oo 

for allT (zM, and all measurable functions F: C([0,oo)) — >■ [0,oo). 
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Proof. It suffices (see e.g. Theorem 14.12 in [21]) to establish (|80f for Fniij) := nr=i MvtJ 
where < . . . < i„ G R and /i, . . . ,/„ e C{,([0,oo), [0, oo)). If i^„(0) > 0, then both sides 
of (80 1 arc infinite. For the rest of the proof, we assume Fn{0) — 0, that is, /i(0) = for at 
least one i e {1, . . . , n}. Wc may even assume fi £ Cc((0, oo)) for at least one i S {1, . . . , n}. 
Otherwise approximate fi monotonically from below with test functions which have compact 
support. In addition, we may without loss of generality assume ti = = T . Otherwise use 
a time translation. 

If Fn vanishes on {rj: rjo = or ryt^ = 0}, then (80) is essentially (78). To see this, 
consider 



(79) 



(77) 



ds Q{dri) 

''-=(F„((rO,>o)ly,„>o)dsQ(d^7) 
(K((i^t)te[o,t„])lyt„>o)'7i(dx). 



E' 



(81) 



Applying (78) with T = tn, reversing the calculation in (81) and substituting s — i„ i— >■ s 
shows that 

/ / Fn{{'nt-s)teii)'^r,-s>o'^m„-s>odsQ{d'q) 

J J — OO 

Fn{{v-t-s)teB)'^v^s>o'iv-t„-s>odsQ{dri). (82) 



We prove (80) with F replaced by Fn by induction on n € IN. The base case n = 1 
follows from a time translation. The induction step n — 1 — >■ n follows directly from (82) if 
/i(0) = = /„(0). If/„(0) >0, then 



K((?yt-s)tgR)lr,_,>0lr,t„_,=0 dsQ{df]) 



/n(0) 
/n(0) 



-F„_i((?7t_s)jg]R)l,,_^>o(l - l,,,^_,>o) dsQ{di]) 



Fn-l{{V-t-s)teB)'^V-s>o{^ 



'-V-t„-s> 



o) dsQidri) 



(83) 



-F'«(('7-t-s)tg]R)lr,_,>ol»7-t„_,=o dsQid-q). 
For the second step we used linearity, applied the induction hypothesis and equation (|82 



and again used hnearity. Adding (82) and (83) proves the induction step in case of /i(0) — 0. 
The remaining case /i (0) > follows from a similar calculation as in ( 83 ) . This completes 
the proof of Lemma [l8| □ 



Lemma 19. Assume [A2l\ and [AE3[ Then 
^'=F{{Yt\^^,^^^)m{dx) - 



F{{7lT-t-s)teiQ,T])ds Qidv) 



(84) 



for all measurable functions F: C([0,T]) — > [0, oo) satisfying F{0) — and all T >0. 
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Proof. Express the speed measure in terms of the excursion measure as in ( 77 1 

E-i^((y,)t6[o,T])™('^2;) 

>oE''-' ( (yOte [o,T] ) 
/ / lr,_,>oF{{r]t-s)te[o.T])'^sQidv) 

J J —oo 

F{{TjT-t-s)te[o.T])Qidv) ds. 



(85) 



— oo 
T 



The last two steps are the Markov property (79 1 and Lemma 18 respectively. 



□ 



With Lemm a |19| in hand, we now reverse time to prove a first version of the Poisson 
approximation (|45|. 



Lemma 20. Assume 



AJ^.l and A2.3. Let c G [0,cx)) be a constant. Let (Yf^g''^)t>s be the 
solution of (44) with Cn{') = c and y,^'*^ = 0. Then 

t f 

ds (86) 



lim NM'fiYf",'') =c f f I{xt^s)Q{dx) 
for all t > and every globally Lipschitz continuous function f on [0, oo) with /(O) = 



Proof. If c = 0, then both sides of (861 are equal to zero. So for the rest of the proof we 
assume c > 0. Let (5 > be arbitrary. Lemma 15 provides us with tojv((0, (5))/A^ ^ as 
N ^ oo. Thus 



lim NM°f{Yfi^^) = lim ..... , 

N^oc -"^ tfl J Ar^oomAr((0,(5)) W->oo7o 



lim / WPf{Y^''^')mN{dx) 



(87) 



=:clim / E^/(y4^'")mjv(da;) +C5 



where Cs is a suitable function of 5. The term Cg converges to zero as (5 according to 
Lemma fT4l 

The speed measure mjv is an invariant (non-probability) measure for {Y^^^'^)t>o, see e.g. 
Theorem V.54.5 in [M]. Thus we may reverse time. As we let N oo, we will exploit that 



(l^^''^)f>o converges weakly to {Yt)t>o. In addition, mi<i{dy) converges vaguely to m{dy) as 
TV — oo as the densities converge. These observations imply that 



hm / EV(r,^'=)m^(d:r) 



N- 



= lim 



<5^0 



f{y)¥y{Y,%''<S)mNidy) 
fiy)Fy{Yt<S)midy) 
f {Xt-s)'i-x-s<sds Q{dx) 
f{xt-s)dsQ{dx). 



(88) 
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The last but one step is Lemma 
theorem together with Lemma 
Lemma [201 



19 



16 



and the last step follows from the dominated convergence 
Putting (87) and (88) together completes the proof of 

□ 



Next we use induction to generalize Lemma 20 to test functions which depend on finitely 
many time coordinates. For this let £t be the following set of bounded functions on 
C([0,T],/) which depend on finitely many coordinates and which are globally Lipschitz 
continuous in every coordinate 



\ WftiVu)- n e N,0 < ti < . . . < t„ < T, /i, ...,/„ are bounded and 



(89) 



globally Lipschitz continuous} 

for every T g [0,cx)). Due to the Lipschitz continuity and boundedness of /i, 
exists a constant Lp ^ (0,oo) such that 



, /„, there 



\F{i^)-F{fi)\<LFY.W.^ 



V77,77ec([o,r],/) 



(90) 



Note that the set £t is closed under multiplication and separates points. Thus the linear 
span of £t is an algebra which separates points. According to Theorem 3.4.5 in [T3] the 
linear span of £t is distribution determining and so is £t- 



Lemma 21. Assume 



All 



and 



N-I are square Lebesgue integrahle and that 
satisfy ( 44 ) . Then 



A2.^ Suppose thatC- [0,T] [0,oo) and that (n- [0,T] 

oo. Let {Y/^;^)t>s 



hm NM'F{iY,'^/),ys) - 

N~¥oo 

for all s > and every function F ^ Et- 



C in Ll^^ as N 
F{{xt-r)t>s)Q{dx) dr 



(91) 



Proof. Let Lp he such that F satisfies (90). Moreover let F be bounded by Cf. Lemma 13 
implies that 



lim iVE°F 



< lim tvlt Ve'^ii; 



hm 7VE' 



3 = 1 



NX 



(92) 



< Lt ■ uCt lim 

N=>-oo 



where Ct is the finite constant from Lemma 
replaced by (. 



\CN{r)-ar)\dr^O 



13 



Therefore, it suffices to prove (91) with ( 



N 



We begin with the case of ^ being a simple function. W.l.o.g. we consider C(i) = 
'^7=i^i^[ti-i,ti)i^) where ci,...,c„ > and tg — s as we may let F depend trivially on 
further time points. The proof of (91 ) is by induction on n. The case n = I has been settled 



in Lemma |20| For the induction step we split up the left-hand side of ( |91[ ) into two terms 
according to whether the process at time ti is essentially zero or not. In order to formalize the 
notion "essentially zero", let (5 > be arbitrary and choose functions ips G C^([0, oo), [0, 1]) 
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such that (f>s jx) — 1 for x > 2S and (j)s (x) = ior x < 6. Furthermore let Ct be the constant 
from Lemma 13 and define F2{ri) :— Jl"=2 fiivu)- 

First wc consider the case that the process is away from at time ti. The fohowing 
estimate shows that we may discard immigration after time ti. The moment estimate of 
Lemma [T3| imphes that 



'ten 



< N]S'>^s(Yt';j)CFnCT-ma.xc,t„ 

iV i^n 



(93) 



> 



for every fixed S > Q. Here we used that K^''' converges to the zero function in distribution 
as — > oo. If the process is essentiaUy zero at time ti (that is 1 — « 1), then 

we may restart the process at time ti in the state 0. The Lipschitz continuity of /i, . . . , /„ 
together with the moment estimate of Lemma 13 provides us with 



UmsupA^E" 



hen 



((^MiOteR) 



< hmsupATE^ (1 - (l)s{Y/;;;;'))nLFCTYf^ 

TV— »-oo L 



uLfCtCi 



N,ci 



(94) 



l-0<5(Xti-r) Xti-r Q{dx)dr. 



The last step is Lemma [20} By the dominated convergence theorem together with Lemma [TB} 
the right-hand side of (|94[) converges to zero as ^ 0. Therefore, using (93), (94) and 
applying the induction hypothesis 



hm nw.'fUy,':,:^) 



AT-i-oo 



lim 7VE° 



lim lim NWP 

<5-i-0JV-s.oo 

+ lim lim E^ 

5-fOAf-i-oo 



lim 

<5-i-0 



J Mxu-r)fi{xu-r)^'''^-''F2[{Yt^t,),^^)Q{dx)dr 
+ /i(0) fcir) I F2[{xt-r\en)Qidx)> 



(95) 



\dr 



C{r) J F({xt-r)ten)Q{dx)dT. 

The last step follows from the pointwise convergence (j)s{x) — )• 1 as (5 — for every x > Q to- 
gether with the dominated convergence theorem and from the Markov property ( 79 ) . Finally 
let C be integrable. Approximate C with simple functions (Cn)nGN- Applying Lemma 13 it is 
straight forward to show that equation (91 ) with C, replaced by C„ converges to equation (91 ) 
as n — oo. □ 



Lemma 22. Assume 



A4--1 and A2.3 Fix s < T. Suppose that CXn- [s,T] — > [0,oo) are 
square integrable and that C, as N oo in the L^-norm. Let iY^g'^ (i))t>s, i < N , he 
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independent solutions of (44 1 satisfying ^/^'''(i) — for every N GlM. Then we have that 

as N ^ oo (96) 



N 



s<t<T 



where 11 is a Poisson point process on [s,T] x U with intensity measure 

Ciu)du(g)Q{d7]). (97) 
Moreover let F be a continuous function on C([s,r],R) satisfying the Lipschitz condition 



\F{rj) - F{fj) \ < LpJ2 



V77,f? e C([s,T],R) 



(98) 



for some s < ti < ■■■ < i„ < some n € N and some Lp G (0,oo). In addition let 
f : I ^ M. be a continuous function satisfying |/(a;)| < Ljx for all x £ I. Then we have that 



lim E 

Af-»-oo 



/ ^ 



i<t<T 



= E 



F 



/(?7t-«)n(dM, dry) 



3<t<T 



(99) 



Proof. Let F C C^((0, oo)) be a countable dense subset of Cc((0, cx))). Tightness of 



N 



i=l 



i<t<T 



TV e N 



(100) 



follows from tightness of 



N 



i<t<T 



N V / € J". 



(101) 



This type of argument has been established in Theorem 2.1 of Roelly-Coppoletta (1986) 
for the weak topology and Cq. Following the proof hereof, one can show the analogous 

argument for the vague topology and C^. Fix / G J" and define 5*^ := Y.^^i /(^m'''(«)) for 
aW s < t < T and TV G N. Note that / is globally Lipschitz continuous. For if G N and 
fixed t G [0,00), global Lipschitz continuity of / implies that 



5f I >if 1 < 



K 



r N 



-E 



< ^ sup ME° 

" MSN 



(102) 



for some constant Lj < 00 and for all G N. The right-hand side is finite according to 



Lemma 



13 



and converges to zero as if — >■ 00. This proves tightness of , G N, for every 
t G [0, 00). For the second part of the Aldous criterion, fix T G [0, 00) and let r^r, G N, be 
stopping times which are uniformly bounded by T. In addition define 



'^-^~x + fiN{x))f\x) + ]^<7l[x)f"{x) Va;G/,tG [s,T],A^GN. 



(103) 

The functions /iat, un and af^ are uniformly globally Lipschitz continuous on the support 
of / according to Assumption A4.1 Therefore there exists a constant C/ G [1, 00) such that 
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\Gff{x)\ <_ Cf{i^+x) and {/(TNfix) < Cjx for all a; e / and all iV € N. For fixed 
T] > and 6 G [0, 1], we use Ito's formula to obtain that 



qN qN \ ^ 



N 



E 



< 3C?E / Cw(^iv 







3C|E 



N i-S 



-NX 



{i)du 



du 



(104) 



< 3C|(5E 



Tiv+5 



3C|(5E 



5 , W 



< (5-3C^ sup 



AfeN Js 



T+1 

[Cm{u)) 



du + (5-6C^ sup E" 



MeN 



M 



for all J < J and for all iV e IN. In the last step we used that (Jp h(u)du) < 6 LJh{u))^du 
for every integrable function h. The right-hand side of ( 104| is finite by Lemma 12 Letting 
^ — >■ 0, the left-hand side of ( 104 1 converges to zero uniformly in € N and S < 5. 

Next we prove convergence of finite-dimensional distributions. Let n e N, / e J^, s < 
ti < ■ ■ ■ < tn and Ai, • • • , A„ > be arbitrary. Using independence we obtain that 



N 



= E 



n[ 

i=l 



Ee" 



(105) 



A^EO 



N 



N 



for all AT e N. Note that 



1 _e- 2:^=1 ^^/(-:') =^(1 



(106) 



for all xi, . . . , a;„ G [0, oo). Applying Lemma 21 to each summand of this telescope sum we 
get that 



lim E 



exp 
E 



C(r)(l - e"^J=i^^^('"^-'))drg(d77)j 



(107) 



3 = 1 
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This proves convergence of finite-dimensional distributions. 

It remains to prove (99), which includes non-bounded test functions. By the previous 
step and by the Skorokhod representation of weak convergence (e.g. Theorem 11.86. 1 in [33] ) 
there exists a version of 

N 



iV G N 



s<t<T 



(108) 



which converges almost surely as — oo. For every K > 2, let ■ [0, oo) — > [0,1] be a 
continuous function which satisfies hxix) — 1 for every x G [;^, and hxix) = for every 
x^[^,2K]. Then /-/ik e Cc((0, oo)). Thus 



N 



t<T 



(109) 



converges almost surely as — > oo for every M G N. The second moment estimate of 
Lemma 



12 



together with |/(a;)| < Ljx for all x G / implies that the argument of F in (109) 
is almost surely uniformly bounded. Consequently taking minimum with M has no effect 
for sufficiently large Af G N and ( 109 ) with M — oo converges almost surely, too. Uniform 
intcgrability of ( 109 ) with M = oo follows from the L^-cstimate 



E 



< 2F{0) + 2i|L|E 



s<t<T 



2n 



N 



3=1 i=l 

N 



(110) 



< 2F 



2L}.Ljn^]S 



sup 

3<t<T 



(E^M^w 



and Lemma 12 Thus ( 109 1 with Af = oo converges also in expectation. The Lipschitz 
condition (98) implies that 

N 



5<t<T 



E 



^ i=l 
n r N 



^((E(/"'^-)(^m'«) 



s<t<T 



3 = i 



N 



<^L^L,E^<f(0( 



< — ^E 
- K 



'-1=1 



sup \YtJ^\ 

s<t<T 



N 

nLpLj sup E[^y,;^^^(^)A- 
s<t<T K 



(111) 



for all G N. Letting N oo and then K — > oo, the right-hand side converges to zero. 
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Thus 



lim E 

Af~>-oo 



= lim lim E 



i<t<T 



lim E 



E 



F 



^((E/(^M») 

^ (E(/"M(^-«) 

u, dri) 



i<t<T 



s<t<T 



(112) 



The last equality follows from the dominated convergence theorem together with Assump- 

□ 



tion 



A2.3 This proves (99). 



4.4 Convergence of the loop-free process 



Recall the loop-free A^-island process from ( 43 ) . The following lemma shows that the loop- 
free TV-island process converges to the virgin island model. 

Lemma 23. Assume lIJl\\A2.3\ \A3.2\ and \A4.S\ Then we have that 

N OO 



J2 J2 V' 



(0 



=1 fc=0 



t<T 



as iV —> OO 



(113) 



t<T 



for every T £ [0, oo). 

Proof. Fix T e [0, cx)). Let C C^((0, oo)) be a countable dense subset of Cc((0, oo)). 
Tightness of the left-hand side of (1131 in G N follows from tightness of 



, W oo 

{(EE/(^"''«)),<,^^^^ 

i=l fc=o - 



V/e 



(114) 



This type of argument has been established in Theorem 2.1 of Roelly-Coppoletta (1986) for 
the weak topology and Cq. Following the proof hereof, one can show the analogous argument 
for the vague topology and Cc- Fix f & and define :— J2iLi SfcLo f{^t''^{^)) 
t <T and A^ S N. For fixed t e [0, oo), global Lipschitz continuity of / implies that 



p(|S'f I >K^< 



K 



-E 



N oo 



i=l fe=0 



< sup E 
A AfeN 



M oo 



EE^^'^w 



i=l fe=0 



(115) 



for some constant Lf < oo and for all AT, A^ e N. The right-hand side is finite according to 
Lemma [t] This proves tightness of ( |114[ ) for every fixed time point. For the second part of 
the Aldous criterion, let tjv, A^ G N, be stopping times which are uniformly bounded by T. 
In addition define H{x) :— Y^^=i SfcLo fi^t) ^^'^ 



N oo ^ N 

g^H{x) := E E (]V E 4"' - + lfc=oMA.(0) + /iAr(x,^))/'(xf ) 



i=l k=0 
N oo 



J = l 



(116) 



EE^^^(^')/"(^') 



1=1 fe=o 
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for all X — {x'^)i<N,ke¥io & l'^ ^ Assumption A4.l| implies that the functions fiN, crjy 

and af^ are uniformly globally Lipschitz continuous on the support of /. Moreover N^pf{0) 
is bounded by 20 uniformly in iV S N. Therefore there exists a constant Ch G [1, oo) be 
such that \g'^H{x)\ < Ch{20 + ^ti T,T=o ^^') ah x = (40^<JV,fcGNo and aU TV e N and 
such that ((/'(TAr)(y)) < ChU for all y G / and e IN. For fixed ry > and S G [0, 1], we 
use Ito's formula to obtain that 



E 



N oo ^TN+S 



(/ g^H{z:!' i-))du + J2i: / {f'■a^){z:!^H^))dBti^) 

2" 



AT 



< sup E 



261 



i=l k=0 

M oo 

sup (2e + Y^Y.z"'{ 

j=l fc=0 



TV oo „5 



(i))(iit 



(117) 



for all 5 < 5 and for all A" G N. The second step follows from (a + 6)^ < 2a^ + 26^ for 

a, 6 S R and from Ito's isometry. In the last step we used that (/^ h{u)duj < /q {h{u)Ydu 
for every integrablc function h. The right-hand side of ( 117[ ) is finite according to Lemma [9] 
Letting (5^0, the left-hand side of (117) converges to zero uniformly in A^ S N and 5 < 5. 



Tightness of (114 1 now follows from the Aldous criterion, see Aldous (1978). 



It remains to identify the limit of the left-hand side of (113) by proving that 



n m N n 

^lim E exp (- E E E /«''(^))) - E exp (- ^ 



i=i 



fc=0 i=l 



y: fivt.^s)) (118) 



A.E 



for all Ai, . . . , A„, for all < ii < • • • < t„ and for all m,n G Nq. Lemma [5] justifies to restrict 
the summation over k to finitely many summands. We prove (118) by induction on m G Nq 



using the Poisson limit (96 1 for independent one-dimensional diffusions. Define F(r]) := 
^jfivtj) for every 77 G C([0,r],R). Note that F satisfies the Lipschitz condition ([25 
for some constant Lp E (0, 00). Let (Yf^'''(i)) ^ , i G N, and (^'j^''' (*)),-^n, « G N, be 



independent solutions of (44) with Y^Q^'ii) 



t>0' ' ■ ' V't,u \-Vt>o' 

Xoil), Y^^'^Hi) = and Cn^) = NfiN{0) for 



every A^ G N and i G N. In addition let (Yt{i))^^^, i G N, be independent solutions of ([s]) 

with Yo{-) = Xoi-). Note that Z.^'°(i) is a solution of (|44| with Z^'°{i) = X^{i). The first 
moment estimate of Lemma [13] implies that 



E 



exp 



N 

E 

1=1 



't>0 



< LpCt„ 



E 



X„^(^f)-Xo(z) 



E 



exp 



E 

i=l 



F 



r N 

E 

1=1 



for all A^ G N where Ct^ is the finite constant of Lemma 13 Letting N ^ 00 the right-hand 
side converges to zero according to Assumption A3. 2 The process (X^''' («)) .^^ in turn 
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is close to except for islands with a significant amount of mass at time zero. 

Formalizing this we use Lemma [13] to obtain that 



lim lim 

X— J-oo voo 



E 



cxp 



N 

E 

= K+1 



t>0 



-E 



exp 



N 

E 

1=1 



t>0 



< lim LpnCt V EfXnfi)] + hm lim E 

i=K+l 



1 — exp 



K 

E 

i=l 



The first summand on the right-hand side is zero ac cordin g to Assumption A3. 2 Note that 



UN M and aN — tr as — oo by Assumption A4.1 Thus (^t^''' (*)) j>o converges in 
distribution to the zero function as N ^ oo for every fixed i € N. Consequently the second 
summand on the right-hand side is zero. Moreover {Yt%'''ii)) converges in distribution to 

(^t(i))j>Q as A^ — > oo for every fixed z S N. These observations imply that 



hm E 

JV->oo 



exp 



N 

E^ 



lim lim E 



exp 



K 

E 

i=l 



(^))*>0 



t>0 



E 



exp 



N 

E 

i=K+l 



hm E 

K^oo 



exp 



K 

E 

1=1 



E 



E 



exp 



i=l 



t>0 



lim E 



exp 



E 

1=1 



(0) 



t>0 



E 



{s,-q)£n<» 



'^^P V E P{{Vt-s)t>a] 

(t,s,r))GV(0) 



The last but one step follows from Lemma 22 with Cjv(0 = A^Mw(O) and C(' 
proves (1181 in the base case m = 0. 



This 



For the induction step m — ^ m -(- 1 note that a version of (Z. 



N,m+1 



on ^r^'"^{j) —■ (Nir), r > 0, is given by the one-dimensional diffusion {Y"q^) 

with vanishing immigration. Thus we may realize {Z^'"'"'^^ {i))t>o by choosing a suitable 



(*))t>o conditioned 



t>o 



version of {Z^'"^{j))t>a, j = 1, . . . , N, and by independently sampling a version of (F/q''") 

whose driving Brownian motion is independent of {{Z^'"'' {j))t>o : j = 1, . . . , N}. Tightness 
together with the induction hypothesis implies that 



i>0 



N m 



E E V'' 

i=l k=0 



(0 



t<T 



E E ' 



as A^ — > oo 



(119) 



t<T 



for every rh < m. Thanks to the Skorokhod representation of weak convergence (e.g. Theorem 
11.86. 1 in [33]), we may assume that the convergence in (1191 holds almost surely. As a 
consequence we obtain that 



N 



t<T 



E 



((.,s,r))eV<") 



as A^ 



oo 



(120) 



t<T 
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holds almost surely. Using arguments from the proof of ( 99 ) one can deduce from this that 

N 



t<T 



(121) 



t<T 



holds almost surely where the total mass of the n-th generation of the virgin island model is 
defined as 

(t,s,r,)eV(") 



(122) 



for every n S Nq. Together with continuity of ('7s)^<y > /q IVs — Vs \ ds, (121 1 implies 
that 

/ |Vzf''"(i)-P;(")|ds^^^O almost surely. (123) 
•^0 ^=l 

That this functional is not bou nded is remedied by a truncation argument. Now the main 
step of the proof is Lemma 22 with CAr(i) := J2iLi ^t^'""!*) ^^.d ({t) :— v/™'' for all t > 0. 
Lemma [22] implies that 

N 



E 



exp 



(-E 



F 



>t<T 



M,m\ 



' AfeN 



= E 



exp(-E^((>^*"0.<.))|(^"'") 

(- / ^^((7/t_,)t<T)n(™+i)(dx,d77))|v(") 



(124) 



exp 



as iV — oo 



almost surely. Here n(™+^) conditioned on V'™-* is a Poisson point process on [0, oo) x U 
with intensity measure 



n('"+i)(ds,d?7)|v('") 



(125) 



= J2 i^s-rds (E) Qidrj) = J2 ^^^''''''^Hds,dT]). 

Due to this decomposition of n.^"^~^^\ we may realize 11^"^+^^ conditioned on V^"*-* as the 
independent superposition of {H^'-'''''^^ : (t, s,?/') S V*-™-*}. In other words, n^™^-'^) is equal in 
distribution to the (m + l)st-generation of the virgin island model. Therefore we get that 

m+l N 



lim E 



N- 



lt<T 



E 



E 



k=0 1=1 

m N 

1™, (-EE [exp (- E F{Z^'"^+\^)) I (Z.^'^-) 

E F({r]t-s)t<T 



E 



fc=0 j=l 

i«^pf"E E ^((^*-)* 

m+l 

exp (- E E Pii'lt-s) 

k=0 (t.s.rjjeVC") 



MeN 



(t,s,r;)eV("'+i) 



which proves ( 118 1 and completes the proof of Lemma 23 



□ 
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4.5 Reducing to the loop-free (A^, /i^v, cr^)-process 

In this subsection, we show that the (A^, (T^)-process with migration levels and the loop- 
free {N, /iAT, cr^)-process are identical in the limit N — > oo. Our proof formalizes the following 
intuition. The individuals of a certain migration level are concentrated on essentially finitely 
many islands. That these finitely many islands are populated by migrants of a different 
migration level has a probability of order As a consequence, all individuals on one fixed 
island have the same migration level in the limit N ^ oo. This intuition is subject of 
Lemma [25] 

First we show that a generation cannot be dispersed uniformly over all islands. To obtain 
this interpretation from the following lemma, assume 

(*) « jf for alH < A^ and some 
time t > 0. Then the cutting operation in (1261 has no effect for A^ large enough. However 
it is clear that the total mass of all individuals with migration level k does not tend to zero 



as A^ ^ oo. Thus A^ 
Lemma 24. Assume 



N,k 



N 



cannot be true. 



A4.L A3. 2 and A4.2 Then any solution of (42) satisfies 



V sup E 



r N 

E 



5^0 



> 



(126) 



for all t G [0,00). The assertion is also true if X^'^{i) is replaced by Z^''^{i). 

Proof. Fix T e [0,oo) and t £ [0,T]. According to Lemma [sj for each e > 0, there exists a 
fco G N such that 



k>kt 



NeK s<T 



N 

sup supE(^Af^'=(z)) <e 



1 



(127) 



Thus it suffices to prove convergence of every summand in (126). In addition, if we forget 
the migration levels in the (A^, /^jv, o'^)-pi'Ocess with migration levels, then we obtain the 
A^- island model. More formally. Lemma [6] shows that 



Af (z) 5] Af '"(0 V*<A^,i>0 



(128) 



m>0 



10 



defines an A^-island model. Recall from (50) and fix AT G N. According to Lemma 
it suffices to prove ( 126 1 with expectation being restricted to the event {f^ < t} for every 
Af G N. Now Lemma 3.3 in 20J shows that, on the event {f^ > t}, (i) is stochastically 
bounded above by y^^"^+^^"(o)_ ^^^^^ Af, AT G N, fc G Nq and (5 > that 



E 



{xri^) 



■ N 

E 

N 

<^E 

4=1 

N 



A (5 1 



^E^o 



4=1 



\ t,o 



>t 



\K+Nhn{0) 



/\6 



■5E 



E ^X^{t)>S 

4=1 



!e 

S 



N 

E- 

i=l 



< A^E' 



K+NfiN{0) 



AS]+Ct 



N 

E 

i=l 



i=l 



- sup E 

6 MSN 



M 



(129) 



[E^o'(*)' 

i=l 
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A3.2 



for some constant C/ < oo. The last step follows from Lemma 13 Next we let N oo and 
(5 — >■ in (129 1. Applying Lemma 22 with (n ■— K + 7V/ijv(0) and using Assumption 
we see that the limit of ( 129 1 as iV — oo and as (5 — >■ is bounded above by 

ft 



\im(K- 



{Xt-r A 5) dr Q{dx) + CtE (Xo(i) 

i=l 



A6 



(130) 



for every if S N and S > 0. The first summand in (|130| is zero by the dominated convergence 
theorem and Lemma 



16 



The last summand converges to zero as (5 — > by the dominated 

convergence theorem and Assumption ^3^2 This completes the proof of (126). The proof 
of ( 126 ) with Xf. ' (i) replaced by '' {{) is similar □ 

Now we prove that all individuals on a fixed island have the same migration level in the 
limit A'' — >■ oo. More precisely, we show that X^'^{i) and X^m^^fc ^t^ '"!*) cannot be big at 
the same time for any i < N or any k g Nq. 



Lemma 25. Assume AJ^-l. A3. 2 and A4.2 Then any solution of (42) satisfies 

N 



N 



:.t<T 



jim sup^E ( J2 E E ) = 



(131) 



for every T e [0,oo). The assertion is also true if X^'''{i) is replaced by Z^''^{i) 



,N.k, 



Proof. We begin with the case s = t. Fix K G IM and recall from (50). By the uniform 
local Lipschitz continuity of fi^, there exists a finite constant Lk such that < Lkx 

for all X < if and all A^ S N. Applying Ito's formula, we see that 



N 



N 



dY.Y.K^k « E <a (^) < 2 E E (^)A^-(o) 



i=l k>0 
N 



m^k 



i=l k>0 



EE<fe«Ef^E^, 



i=l fc>0 
N 



m^k 



i=l k>0 m^k 



N 



\n (j) + LkX • (?) + dMt (i) 



(132) 



E E E <i« ^ E<$ ^0-) + w + d^ri^) 



where {Mi' t>o and {Mi'-^ii)) t>o are suitable martingales for each i < N and k,m e 
Nq. Now take expectations to obtain that 

e(EE<4W E ^ 2m.(0) / E(^^xf4(.)) ds 

^i=l k>0 m^k ' 1=1 fc>0 

+ r^(EE<t««^E E 

"'O ^i=l fc>0 j = lm>0 ' 

+ r ^ (E ^ E 0-) E E <v", (^)) 
+ E (e E <4 « E (^)) 

"'O ^i=l fc>0 m^fk ^ 



(133) 
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for every N G and t < T. Note that the right-hand side is finite. Using Gronwall's 
inequahty, ^Ar(O) < 29/N and J2iLi X]fc>o "'^t^f « (^) ^ K, we conclude that 

-P^(EE<4« E ^ >(4^^ + 2i^^)-e--- (134) 



for every S N. Letting N oo proves ( 131 ) in the case s = t. For the case s < t, apply 
Ito's formula in 



<aV-WE / dX^'"^i^) (135) 

and use similar estimates as above. The case s > t is analogous to the case s < t. □ 

Knowing that there is asymptotically at most one generation on every island, we are now 
in a position to prove that (^X^'''{i)^ and {Z^''' (i)) are close to each other. 

Lemma 26. Assume \l3l\ \M:^ and \l4^ For each iV g N, let 

{(Xf''=(z),Bf(z))^>„:fceNo,z<iV} (136) 



he a solution of (42) and let 

{(^f'(*))t>o: fceNo,z<7V} (137) 



he the unique solution of ( |43[ ) with |(i?|^(i))f>o} replaced by {(_B('^(i))f>o} . Then 

^(EEI^*'^''«-^*'^''«I) ^0 a.7V^^ (138) 



N oo 

^N.,k/,s ryN,k, 
k=0 ' 

for every t G [0, oo). 

Proof. For the time being, assume that a^q and /i^v are uniformly globally Lipschitz contin- 
uous and bounded. The general case will later be handled by a stopping argument. 

Define a;"*" := max(a;,0) and x~ := {—x)^ for all a; € R. We first prove (138) with \x\ 
replaced by x^ and x~ , respectively, separately. As R 9 x i— >■ x^ is not diffcrentiablc, we 
will apply Ito's formula to (f)„ being defined as follows. Let 1 = oq > ai > • • • > a„ > • • • > 
satisfy 

— du=l, / —du = 2,..., / —du = n,... (139) 
1 ^ ia2 ia„ u 

Note that a„ — as n ^ cx). For every n S N, there exists a continuous function : R — >■ 
[0,oo) with support in (a„,a„_i) such that 



2 



nu 



< iJniu) < — for aU M > and / ipn{x) dx = 1. (140) 



Note that tpn{u) < for all u > and n G N. With this, we define the function 



y 



0„(x) := la.>o / dy ip^iu) du, x G R, (141) 
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for every n S N. These functions satisfy 0„ G C^(R), |(/)^(a;)| < 1, (t>'n{x) = ^xya'ipnix), 
4>n{x) < and 0fe(a;)— as A; — ?> oo for every x GM. and n G N. 



Denote the difference process by A* (i) := X, 



N,k/ 



Z^^^{i) for ah i e G, i > 0, fc e No 



and G N. The dependence on N is suppressed for the sake of a more compact notation. 
By definition of (pn, ^ 4>n{x) + a„_i A for ah a: G R and n G N. Thus 



N oo 



1=1 k=0 



M 



, N K ^oo ,Af V 

^ ( E E -^^^ (*)))+ E ^^p ^ ( E («^^- 1 ^ (*)) ) 



(142) 



for aU ii' G N and G N. The first summand on the right-hand side converges to zero as 
— > oo uniformly in G N according to Lemma [H] The last summand on the right-hand 
side converges to zero as TV — oo according to Lemma [24] Consequently 



limsupE( V V fAt*^(i)) ) < limsup limsupE ( V V 0Ar2 (a^ (i)) ) 

^1=1 fe = ' ^4 = 1 fc = ' 

To estimate the right-hand side, we apply Ito's formula to obtain that 
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N K 



(143) 



d^0^.(Af(*)) 



1=1 
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E'^jv^ 



i=l 
N 



\ 



m=0 («) 



N 



(144) 



\ 



m=0 («) 



for every k G Nq. Now we simplify the last summand on the right-hand side using the 
assumption of cr^ begin linear. This linearity implies that x ■ a'jf{y)/y = crff{x) for x = 
X( ' (i) and y = X]m=o"'^t "(O- addition take expectations, estimate |(/)^2(a;)| < 1 and 
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use —(j)jq2{x) ■ X <Q for a: G R. Thus we have that 



, N 
|E[5:0,.(Af(^)) 
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< E 



+ E 



+ E 



[^\l^r\J)\]+^[^\~^N{X^■■'i^))-MZ^■'i^)) 



(145) 



m=0 («) 



for every i > 0, A: e Nq and e N. The last but one summand on the right-hand 



side is estimated as follows. Let 5 > 0. In case of X^'''{i) > 6, apply the inequality 



\-;^fiN{a + b) — ^N{a)\ < 2L^& with a = X^'^{i) and h — X^m^fc ^t" '"'(*)• case of 



X^''^ < S, use |/2Ar(a;)| < Lf^x for all a; > 0. For the last summand on the right-hand side, 
use {y/x~y/y)^ <\x ~ y\ and \ajf{x) — cr^(y)| < Lcr\x — y\ for x,y E [0, oo). Moreover recall 
(x) < jqT^ for all a; > 0. Hence, we get that 
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< e[^ lAr^o-)!] + le[^ |xf - z^-\^\ 
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Xf ■'=(i)><5 



AS 



-E 



for every t > 0, fc e Nq and iV G N. Summing over A: G {0, . . . , K} leads to 

r K N 

EE<^^^(^'«) 



— E 
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r oo TV 



< (l-l-L)E 
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AS 



'-k=0 i=l 



+ iK + l) 



k=0 1=1 
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(147) 



for every t > and iiT, G IN. The last three summands on the right-hand side of (147) 
converge to zero uniformly in t G [0, T] if we first let ^ oo and then (5 — > 0. This follows 
from Lemma 25 and Lemma 24 After inserting (1471 into (143), we see that 

N oo 



limsupEfEEf^'W)^) ^ r(l + i)limsupEfEEk'«)ds 

^,= 1 fc=0 W W^oo V,^^^ ^.^Q I / 



(148) 
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for all t G [0, T]. Note that the right-hand side is finite by Lenima|8] Similarly we obtain ( 148 ) 
with (•)+ replaced by | • |. Finally apply Gronwall's inequality to arrive at (138). 



Next we consider functions cr^ and fipf which are not globally Lipschitz-continuous. For 
each K > choose function cr^ ^ and fJ,N,K which agree with aN and /ijv, respectively, on 
[0, K] and which are uniformly globally Lipschitz continuous and uniformly bounded. Exis- 
tence of such functions follows from the uniform local Lipschitz continuity of aj^ and ^n- Re- 



call the stopping time f^, K,N elN, from (ISO]). The process : i < N,m e Nq} 



agrees with the loop-free (iV, ct|, ^, ^at^k )-process on the event {f^ > t}. Furthermore the 
process {[Xf'"^{i))^^^ : i < N,m € Nq} agrees with an {N, cr'ff x^ /^iv,i<')-process with migra- 
tion levels on {f^ > t}. Therefore 



N oo 

limsupE(^^|xf'=(z) 



1=1 k=0 



for every i^T > by the preceding step. Lemma 
completes the proof of Lemma [26] 

4.6 Proof of Theorem [T] 



10 



(149) 

handles the event {f^ < t}. This 

□ 



Proof of Theorem^ First we prove Theorem [T] imder the additional Assumption |A4.2} This 
will be relaxed later. 

We begin with convergence of finite-dimensional distributions. Recall £t from (89). 
Let F{ri) = YVj=i fjivtj) G £t satisfy the Lipschitz condition (90) with Lipschitz constant 
Lp e (0, oo) and let F be bounded by Cp < oo. Furthermore let the function / : / ^ R 
have compact support in (0, |/|), let / be bounded by C/ < oo and let / be globally Lips- 
chitz continuous with Lipschitz constant Lf € (0, oo). Recall the (A^, /ijy, crjv)-process with 
migration levels from ( [42| . Wc will exploit below that all individuals on one island have the 
same migration level in order to show that 



/ N s , N aa 



1=1 k=0 



t<T 



(150) 



Assuming ( 150 ) we now prove convergence of finite-dimensional distributions. We may re- 



place the (Af, /i^v, (T^)-process with migration levels in ( 150 ) by the loop-free process because 



of Lemma 26 and the Lipschitz continuity of F and /. Hence (150) and Lemma 26 imply 
that 



i=l fc=0 
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, N oo \ , 

j»^<(i:i:/(^."'«))„, =^<( n /( 

\ — 1 ; — n — / ^ I . ^^^■\J 



j=l fe=0 



t<T 



(151) 



The last equality is the convergence of the loop- free (A^, ^at, (T^)-process to the virgin island 
model and has been established in Lemma [23l 

Next we prove (150). According to Lemma [6] if we ignore the migration levels in the 



(A^, /iAT, cr^)-process with migration levels, then we obtain a version of the {N , ,a'j^ j 
process, that is, 



(N ^ , N ca 

i=l - ^ ^ z=l m=0 



t<T 



(152) 
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For proving (1501 we observe that 



fc=0 



< lx"^<S V mSNo + ^ ^ x'^ ^ : 
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fc=0 l^k 



for every sequence (a:'^)fegNo ^ [0: ^-^d every d > 0. Thus we get that 
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(154) 



C(iV,5,t) 



for all G N, t > and i5 > 0. The second summand on the right-hand side converges to 
zero as — > oo according to Lemma[25) The first summand on the right-hand side converges 
to zero as (5 —7- uniformly in G N according to Lemma 24 Using (154) we obtain that 

N 



1 



E/(E 



i—1 m—0 
n r N oo 



t<T 



N 



i=l k=0 



t<T 



^ E^ EE i.-n.)>.|/(E - /K'^xo 



j=l t-i^l k=0 



m=0 
N oo 



Y,c{N,6^t,)+j2^ EEi 

n T- N oo 



x^:''(i)<5 



(155) 



^¥EeEE<''wE<'"« 



n n T- N oo 

+ E ciN, s, t,) + i/ E ^ E E ^ ^ 

i=i j=i ^=1 fc=o 

for all e N and 6 > Q. Letting first N oo and then (5 — > 0, the right-hand side converges 
to zero according to Lemmas 25 and 24 and according to the preceding step. Inserting this 
into ( [1521 ) proves ( [l50| . 

The next step is to prove tightness. This is analogous to the tightness proof in Lemma [22} 
Use the lemmas [T0[ and [T7[ instead of Lemma [12} So we omit this step. 

It remains to prove Theorem |l] in the case when Assumption |A4.2} fails to hold. Fix 
T £ [0,oo ). Let be a bounded continuous function on It follows from Assump- 

tion A3. 2 that \X^\ converges in and thus also in distribution to \Xo\. By the Skorokhod 
representation theorem there exists a version of {X^ {■) : N S N} such that \X^\ converges 
almost surely to |Xo| as A^ — )■ oo. Now the previous step implies that 



AT 



t<T 



N- 



^ E 



H 



E ^v.-. 

(t,M,j;)eV 



t<T 



\Xoi-) 



(156) 
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almost surely. Taking expectations and applying the dominated convergence theorem results 
in 



E 



N 



4=1 ^ t<T 



(t,u,r))eV 



t<T 



almost surely. This finishes the proof of Theorem [TJ 



(157) 

□ 



5 Comparison with the virgin island model 

As in Section |4] we proceed in several steps. First we define the loop-free process. Let 
{(Zf ^(i))t>o: fc e No,i G G} be the solution of 



Sec 



(fc-i) 



{J)m{J,^)~zi''\^)+^,{zi''\^))yt 



(158) 



i eG,k e No, 



where we agree on Z^^^\i) := for i > and i ^ G. We will refer to this process as 
the loop-free (G, m, /i, (T^)-process. The main two steps in our proof of Theorem [2] are as 
follows. Lemma 34 below shows that the total mass of the (G, m, ^, g^)-process is dominated 
by the total mass of the loop-free (G, m, /i, cr^)-process. Lemma 36 then proves that the 



total mass of the loop- free (G, to, /i, cr^)-process is dominated by the total mass of the virgin 
island model. Our proof of Lemma |36] exploits the hierarchical structure of the loop-free 
process. Note that conditioned on migration level k — \, the islands with migration level k 
are independent one-dimensional diffusions. We prepare this in Subsection |5.2| by studying 
the one-dimensional timc-inhomogeneous diffusion 

dY,^f{z) = at) dt - Y^::m -i- (*)) dt + ^J7{ji^)dB^{^ (159) 

where Y^f = x E I and s > 0. The path C G C([0,oo),/) will later represent the mass 
immigrating from lower migration levels. 

The core of the comparison result is the following generator calculation which manifests 
the intuition that separating mass onto different islands increases the total mass. If C = c is 
constant, then a formal generator of {Y^'^]l)t>s is 



g^'^ffix) := {c-x + fi{x))f'{x) + ^a\x)f"{x) xEl 



(160) 



where / G C^(/), see e.g. Section 5.3 in [15 . Recall J-'^}- from (37). 

Lemma 27. Assume \A2.1\ Suppose that c, ci,C2 G / satisfy c < ci + C2. Assume fi to he 
subadditive, that is, ^{x + y) < ^i{x) + /i(y) for all x,y £ I with x + y G I. Let x,y,x + y E I. 
If is superadditive, then 



G^'^f{x + y) < {g^''^f{- + y)){x) + {g^'^'>f{x + -)){y) V/G J-i^lnC^ 



// is subadditive, then 



g^'^fix + y)< (g^^'^fi- + y)) [x) -f [g^^'^fix + •)) (y) V / G n C2 
// is additive, then 

g^'^fix + y)< (g(^^)/(. + y)) (x) + {g^'^^fix + -)){y) V / G n 



(161) 
(162) 
(163) 
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Proof. For / G T^l C^, the first derivative is nonnegative and the second derivative is 
nonpositive. Thus 



g'-''>f{x + y)^ f{x + y)(c-{x + y) + n{x + y)] + ^f'{x + y)a'{x + y) 



1 



<f'{x + y)(ci 



x + ii{x)) + ]^f {x + y)a'^{x) 



(164) 



+ /> + y) (c2 - y + /^(y)) + ^ /" + y)'^'(y) 
^{Q^''^f{- + y)){x) + {Q'^''^f{x + -)){y). 



This is inequahty (161). The proof of inequality (162 1 is analogous. If is additive, then 
a'^{x + y) = cf'^ix) -\-a^{y) and no property of / is needed in the above calculation. 



□ 



Note that the operator on the right-hand side of (|161|) is a formal generator of the superpo- 



sition [Y^y 



Yf^l'^) ,^ of two independent solutions of (|159|). This follows from Theorem 



4.10.1 in Hg. 

We will lift inequality ( 161 1 between formal generators to an inequality between the 
associated semigroups. For this we use the integration by parts formula. For its formula- 
tion, let Gs and Gt be two generators associated with the semigroups {St)t>o and {Tt)t>o, 
respectively. Then, for t € [0,oo), we have that 

rt 

(165) 



Stf-TJ= f Tt_,{Gs~QT)Ssfds 
Jo 



if Ssf e V{Gs) n V{Gt) for all s < t, see p. 367 in Liggett (1985). The idea of using ( [165| ) 
for a comparison is borrowed from Cox et al. (1996). As the generator inequality (161) holds 
for functions in we need to show that the semigroup of {Y^^')t>s preserves J-'^^- This 

is subject of the following subsection. 



5.1 Preservation of convexity 



We write x^ 



(xi 



,Xn) for Xi 



x„ e R and n e N. The i-th unit row vector is 



denoted as for every i e N. RecaU J"!"] from ([3^ 



+ 



(166) 



Lemma 28. For every n € N and f G T^^(jO,oo)) , we have that 

k k 

for all x^ g [0, oo)", all hi, . . . , > 0, all {ii, . . . , i^) G {1, . . . , n}^ and all fc G N. The 
reverse inequality holds if J-'f'J^ is replaced by J-'f'} . 

Proof. The proof is by induction on /c G IN. The base case fc = 1 is trivial. Now assume 
that ( 166) holds for some fc > 1. Fix g:„ G [0, oo)", /ii, /i2 > and (ii, . . . , h-) G {1, . . . , n}'^. 
Applying the induction hypothesis at location -I- X]^=2 ^j^ij ^o the index tuple (zi, ik+i), 
we obtain that 

k k 

/(^„ + X! + ^1^*1 + ^fc+ie^,+i) - /(x„ + ^ hjC,^ + hk+ie,^^^ 



i=2 



i=2 



(167) 



J=2 



3=2 
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The last step is again the induction hypothesis. □ 
Lemma 29. Let n S N, c € [0, oo) and f G J^j'^j^^^ ([0, oo)) . Then the two functions 

) (168) 

and 

f: [0,oo)" ^R,x„ ^ /(a^i,...,x„,c) (169) 

are elements o/ ([0, oo)) . This is also true ij J-j^+ is replaced by J-^ and J-+±, respec- 
tively. 

Proof. The functions / and / are non-decreasing and either bounded or nonnegative. It is 
clear that / is again (i, j)-convex for 1 < i, j < n and that / is (i, ?)-conve x fo r 1 < i, j < n—1. 

at location x„_f_i 



28 



It remains to prove (j, n)-convexity of / for 1 < i < n. Applying Lemma 
to the index tuple (i, n,n + 1), we obtain for all hi,h2 > that 

/(x„ + hid + /l2e„, Xn + /12) = f{^n+l + ^l^i + ^^^Sn + /i2e„+l) \x„ + i=x„ 

> + h2en + /i2e„+i) + f{x„+i + hie,) - /(x„+i)) U„+i=a;„ (170) 

= f{x„ + /i2e„,x„ + /i2) + f{x„ + hie,,x„) - /(x„,x„), 



Lemma 30. A 



ssume 



A2.1 



that is, / is (i, n)-convex. □ 

Let n S N and f € J^^'^i^"'^'' ([0, cxd)) . T/ien t/ie function 

E/(xi,...,x„,y,^r") (171) 

is an element of J^^^(^[0, oo)) for every < s < t. If ^ is concave, then this property still 
holds if J^l^"'^'' is replaced by J-'f'_^^'^ and if J-^^. is replaced by respectively. 

Proof. Fix Q < s <t and n e N. We only prove the case of ^ being concave and / e J^^'!^^'' 
as the remaining cases are similar. According to Lemma [29] it suffices to prove that 

/: [0, 00) ^ R, (xi, . . . , x„+i) ^ E/(a;i, . . . , x„, Y^^f"*') (172) 

is an element of , 00)) . Let {Y^^'^)t>s, a; G /, be solutions of ( |159[ ) with respect 

to the same Brownian motion. It is known that Y^'f < Y^'f^'^ holds almost surely for all 
X < X + h ^ I , see e.g. Theorem V.43.1 in [34] for the time-homogeneous case. Thus the 
function / is again non-decreasing. Moreover / inherits (i, j)-convexity from / for every 
1 ^ *j j !i It remains to show that / is (z, n + l)-convex for 1 < j < n -|- 1. If i < n, then 
(i, n -\- l)-convexity of / at the point (xi, . . . ,Xn, Ytf"^^) shows that 

+ hie., + /i2e„+i) = ^f{x^ + hie„ Y,^f-^' + Y^f-^'+''' - Y^^f^^') 
>W.f{x^ + hie,,Y,^f-^^) +ILf{x,^,Yt^^^^^ (173) 

= f{Xn+l + hiCi) + /(x„+i + /l2e„+i) - f{x^+i) 

for every hi,h2 > and {xi, . . . , Xn) G [0, 00)", that is, {i, n + l)-convexity of / in the case 
i < n. One can establish convexity of 

y^->f{xi,...,Xn,y) (174) 
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as in Lemma 6.1 of [20] (this Lemma 6.1 shows concavity if / is {n + l,n + l)-concave and 
smooth). This step uses concavity of /i. Consequently, / is {n + l,n + l)-convex. This 
completes the proof of / S ([0, oo)) . □ 



Lemma 30 extends Proposition 16 of Cox et al. (1996) pTj. This Proposition 16 is used 
TT] to establish a comparison result between diffusions with different diffusion functions, 
|llj . Using the above Lemma 30 this comparison result can be extended 



see Theorem 1 in 
to more general test functions 



5.2 Decomposition of a one-dimensional diffusion with immigration 
into subfamilies 

Feller's branching diffusion with immigration can be decomposed into independent families 
which originate either from an individual at time zero or from an immigrant, see e.g. Theorem 

1.3 in Li and Shiga (1995). A diffusion does in general not agree with its family decomposition 
if individuals interact with each other, e.g. if the branching rate depends on the population 
size. If the drift function is subadditive and if the branching function is superadditive, 
however, then we get at least a comparison result. In that situation, the diffusion is dominated 
by its family decomposition. More precisely, the total mass increases in the order <f+_ if 
we let all subfamilies evolve independently, see Lemma [33] below. The following lemma is a 
first step in this direction. 



Lemma 31. Assume 



A2.1 



Let x,y, X + y G / and let CiC- [OiOo) 
Lebesgue integrable. If fi is concave and is superadditive, then 



[0,oo) be locally 



< 



t>s 



F + _([0,oo)) 



t>s 



V s > 0. (175) 



// n is concave and is subadditive, then inequality 
If H is subadditive and is additive, then inequality (175 



175) holds with replaced by F_|__|_. 

holds with F_| replaced by . 



Proof Let F„(j/) = fnivt,,-- ■ ,mj e F+_([0,oo)) where /„ e ^|'^([0,^)). We begin with 
the case of C, C being simple functions. W.l.o.g. we consider C,{t) — X^ILi Cil[ti_i,ti)(i) and 
C,{t) = 1]"=!^ Cil[t._j,t.)(t) where ci, c„, ci, . . . , c„ > 0, io = s and t„+i = oo as we may 
let Fn depend trivially on further time points. We will prove by induction on n e N that 



C+C,.x+y\ 



< EF, 



(176) 



For the base case n = 1 additionally assume /i € C^. Approximate a and fi with func- 
tions ai,fii € C°°(R) having the following properties. All derivatives (T;^'°\^p\ fc S N, are 
bounded, /z; is concave and erf is superadditive. Both functions vanish at zero. If |/| < oo. 



then ^ii{\I\) < 



all X e /. Let {Yf^f''')t>s be a solution of (159) with a'^ and /i replaced by af and fii, 



Moreover fJ,i{x) — > fi{x) and <Tf{x) — >■ cr (x) as I 



for 



respectively, and let {Y^^!^^''')t>s be an independent version hereof starting in y G I. Then 
X ^ 5*4/1 (x) ■= F/i(F(^g''^'') is twice continuously differentiable for every t > s, see Theorem 



8.4.3 in Gikhman and Skorokhod (1969). In addition. Lemma 
t G [s,ti]. Consequently, we may apply Lemma 27 



30 



proves Stfi & T^l for all 
to Stfi G Tl'l n C2 for every t G [s,ti] 



and the integration by party formula ( 165 1 yields that 



E 



(177) 
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Now as Z — > oo, (Yf'^^''^''')t>o converges weakly to (Y^'^^'^] 



t>o 



for every x € I, see Lemma 19 in 



Cox et al. (1996) for a sketch of the proof. Therefore letting Z — >■ oo in (1771 proves (176) for 
ri = 1 if /i e C^. The case of general /i € J^i^l follows by approximating /i with smooth 



functions in J-^^l ■ For the induction step n 



1, define 



Note that the induction hypothesis implies that 



V(2/i,...,2/„)G^" 



(178) 



(179) 



and that Lemma 



30 



implies that /„ e J-_ 



(n) 



induction hypothesis, we get that 



< 



,Y 



i.x 



< E 



E 



EF, 



fn+l (■ 



xi+yi,...,x. 



Y 



C.y 



Therefore, using the Markov property and the 

J 

Cv,. \ 11 (180) 



Y^ 



Y, 



t>s 



for all x,y € I satisfying x + y € I. The last step follows from the Markov property and from 
independence of the two processes. This proves (176 1. 

In case of general functions C ^^nd ^, approximate ( and C with simple functions Q and 
0, Z € N, respectively. The process K^g'^ converges in the sense of finite-dimensional distri- 
butions in L^, see Lemma [Tlj and due to tightness also weakly to the process Y.'^f. This 
completes the proof as the remaining cases are analogous. □ 



Lemma 32. Assume [A2l\ and [AEM Then we have that 



N 
i—l 



3<t<T 



n{dy,d'q) 



(181) 



s<t<T 



for all X > and all s < T where H is a Poisson point process on [0,cxi) x U with intensity 
measure Leb (8) Q. 

Proof. The proof is analogous to the proof of Lemma [22} For convergence of finite-dimensio- 
nal distributions use the convergence (131 instead of Lemma 22 Tightness follows from an 
estimate as in (104 1 together with boundedness (see Lemma 9.9 in il9j) of second moments. 

□ 

Finally we prove the main result of this subsection. The following lemma shows that the 
total mass increases if we let all subfamilies evolve independently. In the special case of fi 
and being linear, inequality ( 182 ) is actually an equality according to the classical family 
decomposition of Feller's branching diffusion with immigration. 



Lemma 33. Assume A2.1 Let C- [0,oo) — >■ [0, oo) be locally Lehesgue integrable and let 
X E P If the drift function ^ is concave and the diffusion function is superadditive, then 



t>s 



.([0,oo)) 



T]t-s^idy, drj) 



rjt-uiiidu, drf) 



(182) 



t>s 
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for every s > where H is a Poisson point process on [0, oo) x U with intensity measure 
Leb (8) Q and where 11 is an independent Poisson point process on [0, od) x U with intensity 



measure C(s)ds(8' Q- If H is concave and cr^ is subadditive, then (182 1 holds with F 



replaced by F+-|_. If ji is subadditive and is additive, then ( 182[ ) holds with replaced 

by F+±. 

Proof. The idea is to split the initial mass and the immigrating mass into smaller and smaller 
pieces. Fix s > 0. Let /i be concave and let be superadditive. According to Lemma 31 



t>s 



-([0,oo)) 



0,i 



N 



(183) 



for every N £TN where all processes are independent of each other. Letting TV — >■ cx3 in ( 183 1, 



the right-hand side of (1831 converges to the right-hand side of (182), see Lemma 22 and 
Lemma [32] The remaining cases are analogous. □ 



5.3 The (G, m, /i, (T^)-process is dominated by the 
loop-free (G, m, /i, cr^)-process 

Lemma 34. Assume \A2.1 If ^ is concave and is superadditive, then 



< 



i>0 



F + _([0,oo)) 



OO 



k=0 



{k) 



(184) 



t>o 



If ji is concave and is subadditive, then inequality 

,2 



184) holds with F-| replaced by F_|__|_. 

+- 



If II is subadditive and a is additive, then inequality (184) holds withW^ replaced byW^±. 

Proof. Assume that /i is subadditive and that is superadditive. We follow the proof 



of Lemma 31 and begin with a generator calculation similar to Lemma 27 Let Q-^ and 
denote the formal generators of (G, to, /i, (T^)-process and of the loop-free {G,m, ^i.a'^)- 
process, respectively. Assume that / e ([0, (X))^)r\T^^l ([0, cx))<^) depends only on finitely 
many coordinates. Associated with this test function is 



fc=0 



(185) 



where (a;f ^)jgG,A;eNo S /G'xNo_ N^^g ^hat / e D(^^) = C2(/GxNo), The first partial deriva- 
tives fi :— (g§-)/(X]fcLo ■'''^'^')' * € G, are nonnegative and the second partial derivatives 



fa ■ 



k=0 



I e G, 



(186) 
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are nonpositive. Thus we see that 

oo 



fe=0 



ieG '-jeo fc=o 



(fe) 



oo C30 -1 CXD 



fe=0 



oo p ^ 

^ E E /<E -r-o-' ^) - + M-f ^)) + ^ E 



2/ (/c)^ 



fe=0 



fe=0 SgG 

oo r 



leG 



ieG 



(187) 



E E/<E^ 

k=Q UeG jeG 



m(j, z) 



,(fe) 



for every (x|'°')igG.fceNo G /GxNq^ Now we wish to apply the integration by parts for- 
mula (165). In order to guarantee x i— >■ ]E^/(Xs) e D{Q^), approximate / with smooth 



functions in J^^l, approximate /i and a as in the proof of Lemma 31 and approximate G 



with finite sets. Moreover in order to exploit the generator inequality (187) in the integra- 
tion by parts formula (165), we note that a; E^/(X^) e J^j^l ( sec Lemma 6.1 in [20] ) 
Therefore the integration by parts formula (165) together with inequality (187) implies that 



< E^ 



(fe; 



fc=0 



V/e4'2([0,oo)^ 



(188) 



for all (a;^'^'')igG.feeNn ^ /GxNo_ addition note that {Xt)t>o is stochastically non-decreasing 
in its initial configuration, see Lemma 3.3 in |20j . As stochastic monotonicity is the only 
input into the proof of Lemma 



30 



for all / e J"|'!_+^^([0,oo)'^), we have that 



the assertion of Lemma 30 holds also for {X, 



t)t>o- 



9 (xi, . . . , x„) ^ E-"/(a;i, . . . , x„, Xt) e (/^^ 



Thus, 



(189) 



Using this, the assertion follows as in the proof of Lemma [ST] by induction on the number of 
arguments of F S F^ □ 

5.4 The loop-free process is dominated by the virgin island model 

We show in Lemma [36] below that the total mass of the loop- free process is dominated by the 
total mass of the virgin island model. In the proof of this lemma, we use that the Poisson 
point processes appearing in the definition of the virgin island model preserve convexity in 
a suitable way. This is subject of the following lemma. 

Lemma 35. For every vector z — (zi,...,Zm) G [OjCxd)"*, m G N, let 11^-^ be a Poisson 
point process on a polish space S with intensity measure •^iMi where ^i, . . . , /x™ are fixed 

measures on S. If f € J^^"!*"^^ ([0, oo)) , n € N, then the function 



/: [0,( 



\ n+rn 



R, (3:,z)^E/(x,(5,n(^))) 



(190) 



is an element ofJ^^_ " ([0, oo)) for every measurable test function g : S — >■ [0, c»). Analogous 
results hold if J-^ is replaced by and J-+± , respectively. 
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Proof. The function / is still non-decreasing in the first n variables and (z, j)-convex for all 
i ^ hj ^ n. Furthermore / is non-decreasing in the last m variables as II'^-^ is stochastically 
non-decreasing in z. Fix 1 < i,j < m and /ii,/i2 > 0. Let the Poisson point processes 
n(-\ nC*!^') and n^'s'^j) be independent. Fix z e [0, oo)", x £ [0, oo)" and a measurable 
test function g: S ^[0, oo). Note that H^^' -t- n('*i'=') -I- nC^^-^^) has the same distribution as 
Yii^+hiei+h2ej) ^ Therefore, using (n -I- 1, n -|- l)-convexity of / in the point {x, {g, 11*^-^)), we 
obtain that 



(5, n(^+'^^'='+'^^'=^))) + Mf(^x, (ff,n(^))) 
> E/(x, {g, n(^)) + (.9, nC'^^-))) + Mf(^x, {g, n(^)) + (g, nC'^^^)) 



(191) 



f{x,z + hiei) + .fix,z + h2ej). 



For the last step, note that II'-' + IfC^i'^'^ has the same distribution as n(-+''i*^') and that 
n(^) -|- niMf^j) Jias the same distribution as n^-^^^'^^). This proves {n + i,n + j)-convexity 
of / for all 1 < i, j < m. Similar arguments prove {i,n -I- j)-convexity for 1 < i < n and 
1 <j <m. □ 

Recall the total mass process (v/"'')t>o of the n-th generation of the virgin island model 
from ( 122 I for every n g Nq. 

Lemma 36. Assume \A2.1 and A2.2 If fi is concave and is superadditive, then 



\\\ t \Jk=0,...,koJ t>0 -F + _({0,...,fco},[0,oo)) * )k=Q,...M) 



t>0 



(192) 



for every fcg G Ng. If n is concave and is subadditive, then inequality (192 1 holds with 

F_| replaced 62/F+_|_. If ^ is subadditive and is additive, then inequality (192) holds with 

F_| replaced by F+± . 



Proof. We prove (|192|) by induction on fco G Nq. The base case fco = follows from 

fco 



1^(0) I d__ yio) ^ -y^g apply Lemma 33 for the induction step fco 



1. Fix 



Fn+iix) - /«+i(xi'^\ ■ • • , xit'lV^) e F+_({0, . . . , fco + 1}, [0, oo)) 



(193) 



where fci, . . . , fc^+i G {0, . . . , fco -I- 1} and < si < ■ ■ ■ < s„+i. Let H''^, i G G, C G 
C([0, oo), [0, oo)) , be independent Poisson point processes on [0,oo) x U, independent of 

{Z(°), . . . , ZC^")} and with W-^ = H'^ where has intensity measure 



En'^(du, dry) = ({u)du (g) Q{dri). 



(194) 



Note that conditioned on c'f°\-) := Ejgg the law of z'^^°^^\i) is equal to 

A''o) 1 1 (-C'o) 

the law of F q (*) (defined in ( |159 1) where F q * S G, are independent of each 

other. Thus conditioning on . . . , Z'^'"'"^} and applying Lemma 33 with a; = and s = 0, 
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we obtain that 
E 



k=0,...,k 



' fc=0...feo 



< E 

= E 

< E 



POO 



(d^^'' l).=o....,.v / ^-n^'^'^ ('^-^ ^^)) 



(195) 



The last inequaUty foUows from EisG Cf "^O = EjeG ^- ""^ EieG "^0'' «) ^ l^-""^!, where 
we used X^iGG^^O'*) — ^ from Assumption A2.2 and where we used that Fn+i is non- 



.(feo)| 



decreasing. Next we would like to apply the induction hypothesis. However the right- 
hand side of (1951 depends on {\zj:''°''\)t>o through a continuum of time points and not 
only through finitely many time points. To remedy this, we approximate the Poisson point 
process on the right-hand side of (1951 by approximating {\Z^'^°^\)t>o with simple functions. 
For each N 9 m > max(s„, 2n), choose a discretization {io, • • ■ , ^m^} of [0, m] of maximal 
width ^ such that to = 0, <„j2 = m and {si,...,s,i} C {io, ■ • • , ^m^}- Define h ki if 
ti G {si, . . . , Sn} and k :— otherwise. For a path {xt)t>o G C([0, oo), [0, oo)) , define 



(196) 



Note that (DmX){t) ~^ Xt for every t > as m — > oo. Thus the intensity measure 
jEn^^x**""' ((ju, dry) converges to EII'^ ''° {du,drj) as m — )■ oo. This convergence of the in- 
tensity measures implies weak convergence o f th e Poisson point process H^™-^* to the 
Poisson point process n'^'''°°'. Due to Lemma 35 the function /: [0,00)™^+^ — > R defined 
through 



/(xll«\...,xllf =E/«+i(xl?^\...,X^\y^ J V-uIl''--'''\du,dr,)) (197) 
is an element of J^™_^^([0, 00)). Now we apply the induction hypothesis and obtain that 



E 



< lim E f( ('yf*^'^ 



= E 
= E 



' k=0,...,ko 
/•oo 



?7— «n 



{du, drf) 



(198) 



{iy''')k..^...,.^l j V.^uU-'''\du,d,)) 
T/(*;o+i)A 

1=0,. ...ko' J 



Putting ( 195 1 and ( 198 ) together completes the induction step. 



□ 
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5.5 Proof of the comparison result of Theorem [2] 



Proof of Theorem^ We prove the case of /i being concave and cr^ being superadditive. The 
remaining two cases are analogous. According to Lemma |34] we have that 




^ieA ' *>0 ^iGA fc=0 



for every finite subset A C G. Letting A G, we see that the total mass of the (G, to, /i, cr^)- 
process is dominated by the total mass of the loop-free (G, to, /i, cr^)-process. Now we get 
from Lemma [36] that 



for every /cq € N. Letting — >■ oo, we obtain that the total mass of the loop-free 
(G, TO, /z, (T^)-process is dominated by the total mass of the virgin island model. There- 
fore, the total mass of the (G, to, cr^)-process is dominated by the total mass of the virgin 
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